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ABSTRACT 

We present measurements of the angular diameter distance Da(z) and the Hubble 
parameter H(z) at z = 0.35 using the anisotropy of the baryon acoustic oscillation 
(BAO) signal measured in the galaxy clustering distribution of the Sloan Digital Sky 
Survey (SDSS) Data Release 7 (DR7) Luminous Red Galaxies (LRG) sample. Our 
work is the first to apply density-field reconstruction to an anisotropic analysis of 
the acoustic peak. Reconstruction partially removes the effects of non-linear evolution 
and redshift-space distortions in order to sharpen the acoustic signal. We present the 
theoretical framework behind the anisotropic BAO signal and give a detailed account of 
the fitting model we use to extract this signal from the data. Our method focuses only 
on the acoustic peak anisotropy, rather than the more model-dependent anisotropic 
information from the broadband power. We test the robustness of our analysis methods 
on 160 LasDamas DR7 mock catalogues and find that our models are unbiased at the 
~ 0.2% level in measuring the BAO anisotropy. After reconstruction we measure 
D A (z = 0.35) = 1050 ± 38 Mpc and H(z = 0.35) = 84.4 ± 7.1 km/s/Mpc assuming 
a sound horizon of r s = 152.76 Mpc. Note that these measurements are correlated 
with a correlation coefficient of 0.58. This represents a factor of 1.4 improvement in 
the error on Da relative to the pre- reconstruction case; a factor of 1.3 improvement 
is seen for H. 

Key words: distance scale - cosmological parameters - large-scale structure of uni- 
verse - cosmology: theory, observations 



1 INTRODUCTION 

One of the major challenges in modern cosmology is to un- 
derstand cosmic acceleration, a crucial but perplexing dis- 
covery made through observations of di stant supernovae 
jRiess et al.lll998l ; iPerlmutter et al.lll999l ). This expansion 
is typically attributed to "dark energy" whi ch makes up 
nearl y 75% of the universal energy budget (e.g. iMehta et al.l 
2012). Despite its ubiquitous nature, our understanding of 
its physical properties is poor. It is possible to gain lever- 
age on the dark energy equation of state w, as well as other 
cosmological parameters such as Qk through studying the 
cosmic expansion history. We will use an extension of the 
baryon acoustic oscillations method, namely its anisotropic 
signature, to probe cosmic expansion in this paper. 

Baryon acoustic oscillations (BAO) arise from the 
interactions between matter and radiation in the early 
universe. The story begins with the seeding of primordial 
overdensities composed of dark matter, baryons, photons 



and neutrinos. Prior to recombination, the large density 
of free electrons locks the photons and baryons together 
via Thomson scattering. As the overdensities grow grav- 
itationally, they also heat up and the photons become 
very energetic. The radiation pressure becomes large 
enough to push the photon-baryon fluid outwards in a 
spherical soundwave. The dark matter, which remains at 
the central overdensity, exerts a gravitational restoring 
force on the fluid. The competition between the outward 
push by radiation and the inward pull of gravity gives 



plasma JPeebles & Yul 


197d; Sunvaev & Zeldovichl Il970l; 


Bond & Efstathiou 1984 


; iHoltzmanl 19891; Hu & Sugivamal 


19961; iHu & Whitdl 19961; 


Eisenstein & Hu 1998). The period 



of these waves corresponds to a characteristic spatial scale 
known as the acoustic scale or the sound horizon. This 
scale, ~ 150 comoving Mpc, is the distance traveled by the 
sound wave in the plasma prior to recombination. When re- 
combination occurs, the free electrons are quickly captured 
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which ends the coupling between the photons and baryons. 
The photons stream away, forming the cosmic microwave 
background (CMB), and the baryons are left behind at 
characteristic separations corresponding to the acoustic 
scale. This scale is still measurable in the clustering distri- 
bution of galaxies today , which makes it an ideal standard 
ruler (lEisenstein fc Hul 19981 ; Eisenstein. Hu fc Tegmarkl 
19981; iBlake fc Glazebrookl 120031; lEisensteinl 120031: 
Hu fc Haim anl 120031: iLinderl 120031: ISeo fc EisensteirJ 



20031; iMatsubaral 12004 lAbdalla fc Rawlingsl 120051; 



Angulo et al 



Amendola. Quercellini fc Giallongol 120051; 
2005 : iGlazebrook fc Blakd 120051 : iDolnev. Jain fc Takada 
20061 ). We can use this ruler to measure the distance to 
various redshifts using the clustering distribution of galaxies 
measured from large galaxy surveys. This distance-redshift 
relation depends on the values of the cosmological pa- 
rameters, including those governing dark energy, so it can 
be used to infer the cosmic expansion history and the 
cosmology of the universe. 

To use the BAO method, we must first measure the 
acoustic scale from the clustering of galaxies. This is typi- 
cally done statistically using the 2-point correlation function 
of galaxy separations or its Fourier transform, the power 



on the spherically-averaged statistics dCole et al. 


2005; 


Eisenstein et al.l 20051; Hiitsi 20061; Tegmark et al. 


2006; 


Percival et al 


20071; Kazin et al. 


2010al; Percival et al 


1 201C ; 


Beutler et alj 


201 ll; Blake et al.1 


2011a; Padmanabhan et al. 


20121: Xu et al.ll2012l; Mehta et al.1 20121), which onlv allow 



us to measure the spherically-averaged distance Dy{z) oc 
Da{z) 2 /H(z), where z is the median redshift of the galaxy 
sample. This effectively assumes that the clustering of galax- 
ies is isotropic. Most importantly, the Hubble parameter 
H(z) is degenerate with Da(z) in this measure and hence we 
cannot directly probe the cosmic expansion history encoded 
by H{z). 

The clustering of galaxies, however, is not truly 
isotropic. Anisotropies arise from large-scale redshift-space 
distortions c aused by the line-of-sight velocity of galaxies 
jKaiserll 19871 ) and from assuming the wrong cosmology when 
calculating the 2-point statistics. This second point can be 
used to break the degeneracy between H(z) and Da(z). One 
can imagine that if we assume the wrong cosmology, then 
the BAO will appear at slightly different locations along 
the line-of-sight and transverse directions because the line- 
of-sight distances are measured from redshifts and H(z), 
whereas the transverse distances are measured using the an- 
gular separation and Da(z). 

The use of such anisotropies to measure the 
cosmic expansion history was first proposed by 



lAlcock fc Paczvnskil (| 19791 ). In the case of the BAO, 
the anisotropy can be measured usin g the cluster- 
ing s ignal along different directions ( Okumura et al.l 
20081; iGaztanaga. Cabre fc Hul l201ll; IBlake et al l l2011bl ; 



Chuang fc Wang||201ll ; iKazin. Sanchez fc Blantonll2012l ) or 



by looking at higher order mult i poles of the 2-point statistics 
dPadmanabhan fc White] 120081; iTaruva. Saito fc Nishimichil 



l201ll ; IChuang fe Wanel |2012| ). If the clustering 



isotropic, then all higher order multipoles should be zero. 
However, anisotropies introduce power into the even multi- 
poles (the odd multipoles remain zero due to symmetry). 



This fact can be exploited to infer the values of H(z) and 
D A (z) (see 32J). 

In this paper, we will focus on t he multipole method and 
apply the technique described in IPadmanabhan fc White! 
| 20QSl ) to the 7th data release fDR7; lAbazaiian et al.ll2009l ) 
of the Sloan Digital Sky Survey (SDSS). We present the 
first application of this method to a galaxy redshift survey 
and demonstrate its feasibility. We calibrate our methods 
on 160 LasDamas mocks and perform detailed tests to en- 
sure their robustness. We also apply reconstruction, a tech- 
nique for partially remov ing the effects of non- linear struc- 
ture growth on the BAO |Eisenstein et al.ll2007l ). This tech- 
nique has b een tested extensively on the monopole through 
simu l ations dSeo et a l. 2008; Padmanabhan. White fc Cohnl 



20091; iNoh. White fc Padmanabhanl |2009| ; ISeo et al.1 l201fj| ; 



Mehta et al.l 12011 ) and has recently been applied to 



SPSS DR7 data ( Padmanabhan et al.ll2012l ; IXu et alj|2012l ; 
iMeht a et al. I I2012T T The current work takes a first look at 
how reconstruction affects the anisotropic BAO signal and 
uses this to infer H(z) and Da(z) from the DR7 data. 

In !j2]we present the theoretical background for the mul- 
tipole method. In Sj3]we present our analysis techniques with 
emphasis on the intricacies of the fitting model. In [J4]we in- 
troduce the mock catalogues used and the SDSS DR7 data. 
Sj5]and Sj6] present our fitting results and detailed tests of our 
fitting model and reconstruction technique on the mocks and 
data respectively. We present the cosmological implications 
of our Da{z) and H(z) measurements in iJ7]and conclude in 



2 THEORY 

2.1 Background, Basics and Definitions 

There are two main effects that give rise to anisotropic 
galaxy clustering. The first are redshift-space distortions, 
which arise due to the line-of-sight velocities of galaxies 
such as their peculiar motions within clusters (the Finger-of- 
God effect) o r their coherent infall towards overdense regions 
l)Kaiser]ll987l ). These distort our measurements of cosmolog- 
ical redshifts and hence the line-of-sight separation between 
galaxies. Such effects change the shape of the correlation 
function smoothly with scale (i.e. they have no features and 
are purely broadband in nature). 

Anisotropic clustering can also arise if we assume the 
wrong cosmology when calculating the separations between 
galaxies. Since the distribution of matter is mostly isotropic 
at large scales, artificial anisotropies are introduced by cal- 
culating distances assuming the wrong cosmology as each 
cosmology predicts a unique distance scale. Specifically, we 
calculate line-of-sight separations using redshifts and the 
Hubble parameter H (2) while transverse separations are cal- 
culated using the angular separation between pairs of galax- 
ies and the angular diameter distance Da(z). Both H(z) 
and Da (z) are predicted given a set of cosmological param- 
eters; however, if these do not match the true cosmology 
of the universe, we will measure different clustering signals 
along the line-of-sight and transverse directions. This im- 
plies that the BAO signal in the line-of-sight direction will 
be slightly offset from its location in the transverse direc- 
tion. This is a manifestation of t he Alcock-Paczynski tech- 
nique (|Alcock fc Paczvnskilll97^ ). which uses the measured 
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anisotropy in an object thought to be isotropic to constrain 
the true cosmology of the universe. 

In the past, limited survey volume has made it diffi- 
cult to analyze the differential clustering along the line- 
of-sight and transverse directions. As a result, most BAO 
analyses have been based on the s pherically-averaged (i.e. 
monopole) clustering stati stics (e.g. IPercival et alj|20ld and 
IPadmanabhan et al.l 20121 ). which only allow us to measure 
Dv(z), the spherically-averaged distance to redshift z. This 
quantity is defined as 



D v (z) 



(1 



z) 2 D\{z) 



H(z) 



(1) 



which corresponds to two powers of Da(z) from our trans- 
verse distance measure along the two transverse directions 
on the sky and one power of H(z) from our line-of-sight dis- 
tance measure. However, by measuring the anisotropy, we 
will be able to break this degeneracy between Da(z) and 
H(z), and measure these two quantities separately. 

To measure the anisotropy, we must construct a 
clustering model that includes a parameterization of the 
anisotropic signal. We can then fit this model to the data 
and measure this parameter. Essentially one is presented 
with two choices. The first is to perform fi ts to the trans- 
verse and radial correlation functions (e.g. lOkumura et al.l 
I2008T I and the second is to simultaneously fit the monopole 
and higher order multipoles of the clustering statistics 
jPadmanabhan fc WhitlhoOSf ). 

In this work, we will follow the formalism for the multi- 
pole m ethod based on the work of IPadmanabhan fc White] 
120081 ). We note that the BAO can additionally be shifted 
in an isotropic manner if the assumed cosmology is not the 
true cosmology. Isotropic shifts also occur due to non-linear 
structure growth. We define the isotropic shift in BAO po- 
sition as 

D v (z)/r s 



D v,f(z)/r„j 
D\{z) H f (z) 



D\ f (z) H(z) 



1/3 



(2) 



(3) 



where Dy is defined as above and r s is the sound horizon 
(BAO scale). The / subscripts correspond to a fiducial or 
reference cosmology on which we anchor our measurements: 
all isotropic shifts and anisotropic signals in the BAO are 
measured relative to this fiducial cosmology (see H3.3p . This 
parameterization has been used extensively in past BAO 
studies focusing on the monopole. 

We parameterize the anisotropic BAO signal as e 

[H(z) D A (z) \ ( > 

These parameterizations of a and e are derived from 
isotropic coordinate dilations and anisotropic coordinate 
warpings between the true and fiducial cosmology spaces 
(see Equations 1131 fc 1 14[) . Note that if there is no isotropic 
shift, then a = 1. Similarly, the lack of anisotropy implies 
e = 0. Combining Equations Q & we arrive at 

D A (z) a D A ,f(z) 



H(z)r a 



l + e 
1 



*(l + e) : 



■Hf(z)r aJ . 



(5) 
(6) 



By measuring both the isotropic and anisotropic BAO shifts, 
we can separately constrain the angular diameter distance 
Da(z) and the Hubble parameter H(z) at the median red- 
shift of our galaxy sample z. If we denote the errorbars on a 
and e as a a and a t , and the covariance between them as a at , 
then the errorbars on Da(z) and H(z) can be calculated as 



o~d a h Oh 



This yields, 



9Da 9Da 

dec de 

dH OH 

da de 



dD A 3D 

2 ) ( it ii 

da di 



(7) 



2 

= a- 2 o 2 a + (l + e)- 2 o* - 2a- 1 (1 + e)~ 1 a at (8) 

U A 

2 

= or 2 ol+A{l + e)- 2 o 2 t +4a- 1 (l + e)- 1 a ae (9) 



H 2 
a DAH 
D a H 



= -a- 2 (j 2 t + 2(1 + e)- 2 a 2 - a _1 (l + e) - V a 



(10) 



We will also need a method for distinguishing between 
anisotropies introduced by redshift-space distortions into the 
broadband shape of our clustering statistic and those in- 
troduced through assuming the wrong cosmology (i.e. the 
Alcock-Paczynski signal). There exist simple redshift-space 
distortion models that can be used if we are only inter- 
ested in analyzing the anisotropy in the BAO signal and not 
the details of the redshift-space distortions themselves. Any 
residual inadequate matching between these models and the 
actual broadband shape of the data can be mostly compen- 
sated by includ i ng a few addit ional marginalization terms 
ijSeo et al.ll20Qgl : IXu et al.ll2012l ) such as Equation {51]) de- 
scribed in £13.31 



2.2 Formalism for the Correlation Function 

The clustering of galaxies can be measured using the cor- 
relation function £(r) or the power spectrum P(k). Since 
our analysis will focus on the correlation function, we will 
present the formalism for configuration space here and state 
the analogue for the power spectrum in Fourier spac e, which 
can also be found in IPadmanabhan fc White] ()2008l ). 
We begin with a few basic coordinate definitions, 



2 . 2 
r,i + r ± 



2 2 a 

(i = cos = — I 



(11) 
(12) 



where r is the separation between two galaxies and is the 
angle between a galaxy pair and the line-of-sight direction, 
rii is the separation of the galaxies in the line-of-sight di- 
rection and r± is their transverse separation. In the follow- 
ing, unprimed coordinates will denote the fiducial cosmology 
space and primed coordinates will denote the true cosmology 
space. 

The isotropic dilation (a) and anisotropic warping (e) 
parameters are then defined by 

r\i = a(l + e) 2 r }l (13) 
r± = a(l + e) _1 r_L. (14) 



Substituting Equations (|13[) fc (|14|l into the definitions 
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of r' and fi' as in Equations |TT} & (|12|l . we see that 



«i(1 + eW + (1 + £ ) 



= ar[l + 2eL 2 (ji)] + 0(e 2 ) (15) 

where in the last line we have substituted the second order 
Legendre polynomial Z/ 2 (aO = — l)/2. Also, 

a 2 (l + e) 4 rjj 



/' 



a 2 (l + e) 4 r| + a 2 (l + e)~ 2 r 2 
V 2 + 6e(p 2 - fi 4 ) + 0(e 2 ). 



(16) 

The true 2D correlation function £(r') can be Legendre 
decomposed into multipole moments as 



t(?) = X>(r')MM') 



(17) 



where the L^/(//) are Legendre polynomials of order 
Again, if clustering were perfectly isotropic the £ > mo- 
ments would all be zero. Anisotropy introduces power into 
the even-order multipoles, however the odd-order multipoles 
are always zero due to symmetry. We can substitute Equa- 
tion (|15|l into £p(r') and Equation (|16|l into L(i(fi') and 
write 



&/(cw) + 2eL 2 (/i) 



d&/ (ar) 



Lf/(/i) + 3epi(l - ^ 



dlog(r) 
2 ,dL l i{n) 



(18) 



where we have made a Taylor expansion in both large brack- 
ets. 

Finally, in the fiducial cosmology space, we measure the 
multipole moments 



6(r) = 



2^+1 
2 



g{p)L t {p)dii 



(19) 



+3e 



-£(£- l)(£-2) 



(2£-3)(2^- 1) 
^+1) 



^_ 2 (ar) 



+ 



(2^-l)(2l + 3) 
(<?+l)(£ + 2)(£ + 3) 



£f +2 (ar) 



+2e 



+ 



(2^ + 3)(2^ + 5) 

3£(£-l) <£_ 2 K) 



2(2£- 3)(2f- 1) dlog(r) 
^+1) rf^(ar) 



(2^-1) (2£ + 3) rflog(r) 
3(< + !)(< + 2) <2& +2 (ar) 



2(2£ + 3)(2^ + 5) rflog(r) 



(20) 



where we have used the recursion relation for Legendre poly- 
nomials 



(i + l)L t +i(ji) = (2t+l)iiLt(jM)-iLt-i(jM), 
the derivative relation 

dL t (jJ.) _ (£ + 1)\{j,Li(h)-Li + i(jj,)] 



dfi 



A 1 



and the orthogonality of Legendre polynomials 

-Set'- 



1 2 
L i (fi)L i >(jM)dii= 2£^1 1 - 



(21) 
(22) 

(23) 



Here, 5u< is the delta function. Since we measure the corre- 
lation function using our choice of fiducial cosmology, Equa- 
tion (|20p will form the basis of our fitting template for ex- 
tracting the isotropic dilation (a) and the anisotropic warp- 
ing (e) signal from the data. 

The same r elations can be derived in Four ier space from 
the definitions (|Padmanabhan fc Whltell2008l ) 

Ml = « _1 (l + e)" 2 fc|| (24) 
k' ± = Q" 1 (l + e)fc ± . (25) 

The final equation for the measured multipole moments is 
identical to the configuration space case except f<(r) — > 
Pi{k) and the sign on each occurrence of e is flipped. 

2.3 The Anisotropic Signal 

Although anisotropic BAO information exists in all higher 
order multipoles, its magnitude decreases considerably. 
Hence, for the purposes of this study, we will only focus 
on the monopole (£ — 0) and the quadrupole (1 — 2). The 
monopole and quadrupole we expect to measure according 
to Equation (|20|) are 



2 

£o(r) = (,o(ar) + -e 



3£ 2 (ar) + 



d&(ar) 



dlog(r) 



(26) 



&(r) = 2e 



ri£o(ar) 
dlog(r) 
4 

V 



+ 



, 6 V , , 4 d£ 2 (or) 
l + =c)fc» («r) + =e ^ \ 
7 / 7 dlog(r) 



5^4 (ar) + 



d£ 4 (ar) 



(27) 



dlog(r) 

and will form the basis of our analysis. Here, £4(7*) is the 
hexadecapole (£ = 4). 

Figure |l(a)| shows variations of the expected monopole 
(Equation I26|) . the transverse and radial correlation func- 
tions, and the quadrupole (Equation I27p with e for a linear 
theory based model. We have included large-scale redshift- 
space distortions (the Kaiser effect) so that the quadrupole 
becomes non-zero, but not the Finger-of-God (FoG) effect. 
In Fourier space, this model is given by Equation (|28[) with- 
out the F(k, /i, £ s ) term. The transverse and radial corre- 
lation functions were calculated as £0 + £2(^)^2 for /i = 
and 1 respectively. Note that taking the difference between 
these yields the quadrupole. 

One can see that the sensitivity of the monopole to e is 
quite low. However e does cause shifts in the BAO position 
in the line-of-sight and transverse correlation functions. The 
line-of-sight shift is larger than and opposite in direction to 
the transverse shift. This causes the BAO feature in the 
quadrupole at ~ 110/i -1 Mpc to move with e, in addition to 
the shift caused by the isotropic dilation a. Hence we see 
that the quadrupole can be used to obtain a measurement 
of the anisotropic BAO signal via e. 

2.4 A Non-linear Model 



In Figure 1(a) discussed in the previous section, we assumed 



linear theory including the Kaiser effect for £o(r) and £ 2 (r). 
However, in order to model actual observations with fidelity, 
we must also account for the FoG effect and non-linear struc- 
ture growth. This section details a plausible model that in- 
cludes all of these effects and will be used in our fitting 
procedure described in ^3.31 to measure a and e. 
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(a) The monopole (left) and quadrupolc (right) we expect to measure in the presense of anisotropic clustering 
according to Equations l|26|l & l|27|l for a linear theory based model including the Kaiser effect. The Kaiser effect 
gives rise to the BAO bump near 110/i -1 Mpc in the quadrupole. One can see that the monopole is insensitive to 
e. However, e works to change the position of the line-of-sight and transverse BAO features in opposing directions 
with the line-of-sight having a more prominent shift. The differential nature of these shifts moves the quadrupolc 
BAO, compounding with any isotropic shifts. Hence, we expect the quadrupole to be sensitive to an anisotropic 
BAO signal. 
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(b) Variation of monopole (left) and quadrupolc (right) models including the Kaiser effect and a full non-linear 
treatment of FoG and anisotropic E n i . The solid black line in this and the similar plots following always corresponds 
to the fiducial model parameters Ex = 6/i _1 Mpc, Ey = 10h~ Mpc and E s = 4/i Mpc with /3 = 0.35 (center of 
the /3 prior in our fits). The monopole is again affected very little by e. The fiducial quadrupole model picks up a 
crest-trough-crest structure at the BAO scale due to the differential broadening of the line-of-sight and transverse 
BAO signals by Kaiser, FoG and anisotropic Ejj. We again see that the anisotropic warping parameterized by e 
works to shift the location of the quadrupole BAO. In addition, it can adjust the relative amplitude of the crests, 
e is the only parameter that can shift the BAO in the quadrupole while leaving the monopole BAO unaffected; the 
isotropic shift a changes the BAO position equally in both. 

Figure 1. Variation of our models with e for a linear theory based model including the Kaiser redshift-space distortion (a) and a full 
non-linear model including FoG and Kaiser redshift-space distortions as well as anisotropic E n i (b). e parameterizes the amount of 
Alcock-Paczynski anisotropy, which, if there was none, would be equal to 0. The left panel shows the monopole (black), the transverse 
correlation function (red) and the radial correlation function (blue), where the difference between these latter two yields a measurement 
of the quadrupole. Solid, dashed and dotted lines arc defined as in the plot legend of the right panel which shows the quadrupole. Note 
that the quadrupolc BAO feature is much weaker in the more realistic non-linear model. 



In Fourier space we can write the following template for 
the 2D non-linear power spectrum 

P t (k,jJ,) = (l + P^) 2 F(k,n,T, 3 )P dv (k,n) (28) 



where 



(29) 



(1 + fcVE!) 2 

corresponds to a streaming model for the FoG effect and the 
(1 + Pfi 2 ) 2 term corresponds to the Kaiser model for large- 
scale redshift-space distortions. Here E s is the streaming 
scale and is typically ~ 3 — 4/i _1 Mpc. 



The de- wiggled power spectrum Pdw(fc,A t ) is defined as 



Pdw(fc,M) = [flm (AO - Pnw(fe)] 

fcV 2 Ej{+fc 2 (l-ti 2 )E 2 



• exp 



~t~ -fnw (&) 

(30) 



where Pn n (k) is the linear theory power spectrum and 
Pnw(fc) is a power spec trum without an acoustic signature 
|Eisenstein fc Hull 19981) . £„ and Ex are the line-of-sight and 
transverse components of E n i, i.e. E 2 ! = (Em-|-Ex)/2, where 
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E n i is the standard Gaussian damping of the BAO used 
to model the de gradation of the signal due to non-linear 
structure growth ( Eisenstein. Seo fc White~ll2007l ). Here, the 
damping is anisotropic due to the Kaiser effect. 
The multipoles of this template are then 



21- 



Pt(k, n)Li(fi)dfi, 



(31) 



which can be transformed to configuration space using 
fc 3 dlog(fc) 



6,t(r) 



2tt 2 



-Pi,t{k)ji{kr) 



(32) 



Figure 1(b) shows the variation of our non-linear 



monopole and quadrupole models with e while Figure [2] 
shows the variations with other parameters. The variation 
with a is not shown as its role is well understood: a works to 
shift the BAO feature around equally in both the monopole 



and quadrupole. Figure 2(a) shows the variations with Ej 



an d Ej . Figure 2(b) shows the variations with E s and Fig- 
ure 2(c) shows the variations with /3. In these plots, the solid 



black line always corresponds to the fiducial model parame- 
ters E± = 6/i _1 Mpc, E|| = 10/i _1 Mpc and E s = 47i" 1 Mpc. 
P is set to the center of the prior, 0.35, unless indicated oth- 
erwise. These fiducial parameters will be discussed in more 
detail in £13.31 

We see that the monopole model is only weakly af- 
fected by varying these parameters. E n i has an immedi- 
ately obvious effect on the BAO peak since it is designed to 
damp the BAO to model non-linear evolution. The change 
in the peak is only significant with a large change in E n i. 
The Ex = E|| = 8h~ Mpc case and the fiducial case have 
very similar E n i values. Hence, we see little difference be- 
tween the monopole models in these two cases. However, 
the Ex = 4/i _1 Mpc, E y = 7/i _1 Mpc case corresponds to 
E n i ~ 6/i -1 Mpc which affords a weaker smearing of the 
peak. Large E 3 values can also weakly damp the monopole 
BAO, causing slight modifications to its shape. /3 appears 
to have little effect on the monople and e has no effect. The 
only parameter that can shift the monopole BAO position 
is a, which also shifts the quadrupole BAO equally. 

Looking at the quadrupole model in the fiducial case, 
we see that the BAO (at r ~ 110/i -1 Mpc) looks different 



from the case presented in Figure 1(a) There we saw a sin 



gle bump at the acoustic scale due to the Kaiser anisotropy 
(1 + /?/i 2 ) 2 . Including FoG and anisotropic E n i introduces 
more structure at the BAO scale. We see that in the fiducial 
model, a dip has appeared at the center of the linear-theory 
peak, creating a crest-trough-crest structure. This is the re- 
sult of anisotropic E n i, E s and P differentially broadening 
the BAO peak in the radial and line-of-sight directions. We 
see that the radial BAO is wider than the transverse BAO, 
but the radial BAO peak has more contrast. Subtracting the 
two therefore yields the observed crest-trough-crest shape of 
the quadrupole near the BAO. 



Figure 2(a) shows that changing Ex and Em gives rise 



to crests and troughs near the BAO scale in the quadrupole. 
Taking E n i to be isotropic and non-zero completely removes 
the trough, leaving a single peak that is broader than the 
linear-theory case. Changing the values of Ex and Em (but 
keeping their ratios roughly the same) alters the structure 
of the peaks and changes the crest-trough contrast. 

Changing E s has the most prominent effects at small 



scales as expected, since the FoG is strongest there. It also 
causes a noticeable change in the quadrupole BAO signal 
which is partially degenerate with the effects of anisotropic 
E n i, i.e. it can also adjust the crest-trough contrast and can 
eliminate the trough entirely (E 3 = 0ft _1 Mpc). Leverage on 
this parameter mostly comes from small scales, where E s 
has a significant effect on the quadrupole shape. 

Changing p shifts the overall magnitude of the 
quadrupole at large scales. Hence, we see that varying Ex, 
E||, E s and p only affects the shape of the quadrupole, not 
the BAO position. 

Only e and a can shift the location of the BAO. While 
a shifts the BAO position equally along all directions, Fig- 
ure 1(b) shows that e shifts the radial BAO position more 



than and in the opposite direction to the transverse BAO 
position. Therefore changing e will cause the quadrupole 
BAO position to change in addition to the shift induced 



by a. We see from 1(b) that e can also adjust the BAO 
shape so it is not completely non-degenerate with the other 
parameters. We emphasize that Ex, Ey, E s and ft cannot 
shift the quadrupole BAO, they merely work to change the 
BAO shape, e is the only parameter that can change the 
quadrupole BAO position without changing the monopole 
BAO position, so its effects should be detectable. 

The above observations of the model quadrupole be- 
haviour can be re-cast into Figure(3] The panels of this figure 
show the derivatives of the quadrupole model with respect 
to each parameter (Ex - top left, Ey - top right, E s - mid- 
dle left, P - middle right, a - bottom left and e - bottom 
right). The dashed line marks the acoustic scale. The plot- 
ted derivatives show the variation of the model with these 
parameters and are especially interesting near the acoustic 
scale. 

One can see that the behaviour of the Ex derivative 
near the BAO scale is exactly opposite to the Em case. In the 
former we see a down-up-down structure and in the latter we 
see an up-down-up structure. The E s derivative shows sim- 
ilar behaviour to the Ey derivative near the acoustic scale. 
However, its small-scale behaviour is much different. The p 
derivative shows only a single peak near the acoustic scale. 
Note that all of these derivatives are symmetric about the 
acoustic scale. 

The a and e derivatives are different from the others in 
that they are both anti-symmetric about the acoustic scale. 
This reflects their ability to shift the BAO feature. The dif- 
ference between a and e lies mainly in the monopole where 
a can significantly shift the BAO but e cannot; e only shifts 
the BAO in the quadrupole. Near the acoustic scale, the 
quadrupole a derivative shows a large amount of structure 
while the e derivative shows a simple up-down structure. 
The crests and troughs of both of these will be partially de- 
generate with the other parameters despite their opposite 
symmetries near the acoustic scale. However, if we place 
well-informed priors on the other parameters, we limit the 
model from exploring these degeneracies thereby obtaining 
reasonably robust measurements of e. 

A similar conclusion can be reached by noting that the 
first three cases look like the derivative of a Gaussian with 
respect to its width. The P case looks like the derivative of 
a Gaussian with respect to its height and the e case looks 
like the derivative of a Gaussian with respect to its center. 
All of these behaviours are different. 
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(a) Variation of the monopole (left) and quadrupolc (right) models with Ex and Em. The BAO peak in the monopole 
can be affected by E n i as expected since this parameter is used to model the smearing of the BAO due to non-linear 
structure growth. In these plots we have demonstrated the effects of going to an isotropic E n i of roughly the same 
magnitude as the fiducial case (dotted line) and a smaller E n i (dashed line). We see that going to an isotropic E n i 
has little effect on the monopole, however, it completely eliminates the trough feature at the BAO scale in the 
quadrupolc. Going to a smaller E n i makes the peak appear sharper in the monopole as expected. In the quadrupolc, 
it alters the structure of the peaks and reduces the crest-trough contrast. 
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(b) Variation of the monopole (left) and quadrupole (right) models with E s . The BAO feature in the monopole can 
be slightly broadened by a large E s . In the quadrupole, the effects of E s are partially degenerate with E n i in that 
it can alter the crest-trough contrast and can also completely eliminate the trough (E s = 0h~ J Mpc). However, the 
effects of E s are much stronger at small scales due to its large influence on the radial and transverse correlation 
functions, giving us leverage on this parameter. These variations are not surprising since the FoG effect is most 
pronounced at smaller scales. 
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(c) Variation of the monopole (left) and quadrupole (right) models with /3. Again we see that the monopole is 
not significantly affected by /3. The line-of-sight and transverse correlation functions experience large changes in 
amplitude with /3, leading to amplitude differences in the quadrupole at large scales. 



Figure 2. Variation of the non-linear monopole and quadrupole models with different model parameters: E n i (a), E s (b) and (c). 
Comparing the behaviour of these parameters to e (Figure [l(b)[ l indicates that the various model parameters have mostly different effects 
on the quadrupole. While all of these parameters can affect the shape of the quadrupole, only e can change the quadrupole BAO position 
separately from the monopole BAO (a changes both in lock-step). Hence we expect that the effects of e should be detectable. 
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Figure 3. Derivatives of the model quadrupole with respect to T,±_ (top left), Sy (top right), S s (middle left), /3 (middle right), a 
(bottom left) and e (bottom right). The plotted derivatives illustrate how the model changes with these various parameters and is 
especially interesting near the BAO scale marked by the dashed line. Note that near the acoustic scale, the S^, E|| and S s cases look 
like derivatives of a Gaussian with respect to its width. The /3 case looks like the derivative of a Gaussian with respect to its height. 
The e case looks like the derivative of a Gaussian with respect to its center. These behaviours are all different. We see that the Ex and 
Eii derivatives are similar in nature at the acoustic scale but opposite in sign. The Ey and E s derivatives, however, are of the same 
sign and show the same up-down-up structure near the BAO scale, but differ at small scales. The /3 derivative only shows a single peak 
near the acoustic scale. We also see that the Ex, ^||, an d P derivatives are symmetric about the acoustic scale while the a and e 
derivatives are anti-symmetric. The a derivative has the most structure near the acoustic scale. The e derivative shows a simple up-down 
structure. Despite their opposite symmetries near the acoustic scale, the various crests and troughs of the a and e derivatives will be 
partially degenerate with the other parameters. However, given reasonable priors on the other parameters, the model will not be allowed 
to explore these degeneracies and we will recover robust measurements of e. 



Measuring Da and H using BAO 9 



2.5 Covariance Matrix Formalism 

Assuming that the overdensities are drawn from a Gaus- 
sian random field, the covariance between two multipole mo- 
ments £ and £' is 



f K sin(x') , 
si(x) = / V^-dx 



and 



2{2£+\){2l' + f) 



jo(fcr) = 



sin(fcr) 
kr 



(41) 



(42) 



7 



V 

fc 3 dlog(fc) 
2tt 2 



(33) 



Note that due to symmetry, we have Cij(^o(fi), £,2(rj)) = 
Cij (^2(ri), £o(rj)). Finally, the last case has the form 



50 f fc 3 dlog(fc) 



V 



where V is the survey volume, the je(kr) are the spherical 
Bessel functions of order £ and Pf e ,(k) is denned as 



2tt 2 



Pu>{k) 



P(fc, At ) + 



L c (fj,)L e ,(fi)dfi. (34) CyC&W.&'fa)) = 



Collectively, we can then write 

2(2£ + 1)(2£' + 1) 



j2(kr t )J 2 (kr 3 )Pi 2 (k). 

(43) 



Here, P(k, fj,) is the 2D power spectrum and n is the mean 
number density of galaxies. The 1/n term corresponds to 
Poisson shot-noise. We ignore the redshift dependence of n 
above, however, we will describe a method that allows its 
inclusion in £13.31 

In practice, the correlation functions we measure are 
binned and hence, we must also account for this binning in 
the Gaussian covariance matrix. If we calculate the correla- 
tion function in bins with lower bounds r\ and upper bounds 
r 2 , the binned covariance matrix is 



Cy(&(ri),&/fa)) = 



2(2£ + l)(2g' + l) 



V 



2 , dfl 
r ar 



I 



4tt 

fc 3 dlog(fc) 



.3, 



1 3 7 - 3 7-. 3 

i2 'il ' j2 

, 2 

r ar 

47T 



2tt 2 



j e (kr)j e ,(kr')P? e ,(k). 

(35) 



In the case of the monopole and quadrupole relevant to this 
study, the above expression can be split into the £ — £' — 
case, the £ — and £' — 2 (or vice versa) case and the 
£ = £' = 2. T he first case was shown to have the form 
jXu et alj|20l3 ) 



C;j(£o(r;),f (rj)) 



V 



fc 3 dlog(fc) 



2tt 2 



where 



and 



Jo(kr) = 



ji(kr) = 



sin(fcr) 
(kr) 2 



Mkn)J {krj)P% {k) 
(36) 

(37) 



2 ji(kr) 



i(kr) 



kr 



(38) 



Here, [/(x)]^ is standard notation for /(a) — f(b) for any 
function /. The second case has the form 

CyttoC-iUafa)) = / k * d ^f k) Mkn)J2{k r] )PUk) 



V 



where 
Mkr) 



3si(kr) 1 / 3r . . 2 . N 

~~ k 3 k\T J °( ) + r Jl( - 1 



(39) 



(40) 



V 

fc 3 rflog(fc) 
2tt 2 



Mkr^Jvikr^Pfeik). 

(44) 



These forms for the binned Gaussian covariance ma- 
trix form the basis for deriving a covariance matrix for our 
data. We extend the method for approximating the mock 
covariances using a modified for m of the binned Gaussian 
covariance matrix as described in IXu et al.l ||2012] ). This will 
be outlined in i]3.3l 



3 ANALYSIS 

3.1 Reconstruction 



Non-linear structure growth degrades and shifts the acoustic 
peak. Reconstructi on was initially propo sed to partially re- 
move these effects |Eisenstein et al]|2007l '). In the linear the- 
ory description of structure growth, overdensities are small 
and hence remain largely in place as they accrete more 
matter and grow. However, at low redshifts, this descrip- 
tion becomes increasingly less suitable as some overdensities 
grow to masses large enough that they begin exerting sig- 
nificant gravitational pulls on each other. This gives rise to 
pairwise relative velocities between particles separated by 
~ 100/i _1 Mpc. These coherent flows that form large-scale 
structure are the dominant source of smearing of the BAO 
signal. The peculiar motions of particles within a gravita- 
tionally bound structure are subdominant. Reconstruction 
attempts to undo these coherent motions in the matter den- 
sity field and arrive back at something that more closely 
resembles linear theory. This translates into a sharpening 
up of the acoustic peak in the correlation function which al- 
lows us to gain a better centroiding of its location and hence 
a more precise measure of the acoustic scale. It also removes 
some of the shifting of the BAO due to non-linear structure 

;rowth, which has been shown to be a ~ 0.5% effect at 2 = 

Seo et al.|[2010l ; iMehta et al.ll201ll ). 

Our reconstruct i on alg orithm is described in detail in 
IPadmanabhan et alj (|2012l ). A more pedagogical discussion 
of reconstruction can also be found there. This method is 
extended from the origina l recon struction algorithm pro- 
posed in lEisenstein et al. ll2007h and its theoreti c al ba - 
sis is established in IPadm anabha n. White fc Cohnl (|2009h 
and iNoh. White fc Padmanabhanl (|2009h . A basic outline is 
given below. 
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The ultimate goal of reconstruction is to infer the mat- 
ter displacement field that arises due to non-linear struc- 
ture growth from the observed galaxy density field. Then, 
we can shift the galaxies back along their inferred displace- 
ment vectors to place them where they would have been 
in linear theory. This is simple if we consider only the first 
order displacements "J/. In this case 



V 



* + /JV ■(*,§) = -% 
o 



(45) 



where $ s = ■ s is the displa cement in the line-of-sight 
direction l|Nusser fc Davislll994l ). 5 ga i is the galaxy density 
field, b is the large-scale galaxy bias (which is roughly con- 
stant) and hence 5 ga i/fo is an approximation of the matter 
density field. The second term in this equation arises due 
to large-scale redshift-space distortions caused by the co- 
heren t infall of galaxies towards overdense regions (|Kaiserl 
1 19871 ). This implies that as an additional bonus, reconstruc- 
tion can also correct for linear redshift-space distortions. 
The j3 parameter governs the amount of anisotropy intro- 
duced by the Kaiser effect. It is defined as /3 = f/b where 
/ = dlog D(a)/d log(a) ~ fl m (a) 0,55 is the linear growth 
rate, D(a) is the linear growth function and a is the scale 
factor (|Caroll. Press fc Turne3ll992l : lLindeill2005l ). If we as- 
sume that Vl" is irrotational (i.e. curl-free), we can write 
^ = V0. After selecting appropriate values of b and /3, we 
can solve for the scalar field <j) an< A then the displacement 
field "J/ using a finite difference approach. 

Areas that do not fall within the survey or are masked 
out by the survey need to be accounted for as the gravita- 
tional potential (and hence the displacement field) is sen- 
sitive to these regions. Our algorithm fills in these exter- 
nal regions by embedding our survey into a volume that is 
sufficiently larger than the survey volume before the dis- 
placement field is calculated (we use 200/2,~ 1 Mpc padding). 
It then smooths over the density field using a Gaussian 
with a 15/i -1 Mpc smoothing length and generates a con- 
strained realization o f a Gaussian density field that matches 
the o bserved density (|Hoffman fc R ibak 199l1; IZaroubi et alj 
Il995l ). This fills in the masked and unobserved regions. 



3.2 Computation 

The computation of our correlation functions is tied to our 
reconstruction algorithm. We bin our correlation functions 
in 3h~ Mpc bins starting at 2.5h~ Mpc and going up to 
197.5/i _1 Mpc. A list of the steps involved is given below. 

I) Obtain a set of randomly distributed points that have 
the same angular and radial selection function as the survey. 

II) Compute the unreconstructed correlation func- 
tion from the data using the Landy-Szalay estimator 
|Landv fc Szalavlll993f ). 



DD(r, fi) - 2DR(r, y) + RR(r, fi) 
RR{r, fi) 



(46) 



where DD, DR and RR are the number of galaxy- 
galaxy, galaxy-random and random-random pairs that 
are separated by r and ji. We apply FKP weighting 
l|Feldman. Kaiser fc Peacockll 19941 ) for each object as 



(47) 



where n(zi) is the number density at the redshift of the 
object Zi and P(ko) = 40000/i -3 Mpc 3 is the approximate 
value of the power spectrum at the BAO scale. 

III) Estimate the galaxy bias 6 and the anisotropy pa- 
rameter /3 from the unreconstructed correlation function. We 
use fiducial values of b = 2.2 and (3 = 0.3. 

IV) Embed the survey in a larger volume and smooth 
the density field using a Gaussian (again, we use a smoothing 
length of 15/i~ 1 Mpc). Generate a constrained Gaussian re- 
alization matching the observed density to fill in the masked 
and unobserved regions. 

V) Estimate the displacement field Vl" using Equation 
(|45|l and shift the galaxies by -$-/($ s s) to partially undo 
the effects of non-linear structure growth and large-scale 
redshift-space distortions. This is the essence of reconstruc- 
tion. 

VI) Obtain another set of randomly generated particles 
with the same radial and angular selection function as the 
survey. Shift these by — \E" and denote as S. 

VII) Compute the reconstructed correlation function 
using the Landy-Szalay estimator 



£(r,M) = 



DD(r, fi) - 2DS(r, fx) + SS(r, fi) 
RR(r,n) 



(48) 



3.3 Fitting 

We construct models of the monopole and quadrupole in a 
fiducial cosmology to measure the position of the BAO in 
the data relative to the model (parameterized by a) and the 
degree to which it is anisotropic in the data (parameterized 
by e). We base our fitting templates for the monopole and 
quadrupole on Equations ()26|) & (|27[) . That is, we define 
fitting models of the form 



£o,m(r) 
6>,m(r) 

where 



C 2 C , v . 2 I" d£ 2 ,t( ar ) 
B £ ,t{ ar ) + 7 e 3^2, t(ar + ' 

alog(r) 



+ A (r) 



(49) 

+ U^M + ^^\+Mr) (50) 
7 L olog(r) 



Mr) 



at,i 



+ 



+ ae,3- 



(51) 



l + n(zi)P(k ) 



The Ai(r) are composed of linear nuisance terms used to 
marginalize out broadband effects such a s scale-dependen t 
bias and redshift-space distortions as in IXu et al.l (|2012l) . 
The Bo term adjusts the amplitude of the monopole tem- 
plate £o,t- We infer the galaxy bias b from the multiplicative 
offset, b 2 , between this template and the measured correla- 
tion function at r = 50/i~ 1 Mpc. We then use this offset to 
normalize the full models, £o,m and (,2,m in Equations (|49l) & 
(|50|l . to the data. This ensures that Bq ~ 1 as it is primar- 
ily the monopole fit that sets this term and e is very small. 
In practice we perform our fits in the non-linear parame- 
ter log(Bo) to prevent B 2 from going negative which would 
be unphysical. We adopt a Gaussian prior on log(£?o) with 
standard deviation 0.4 and centered at to prevent B 2 fro m 
wandering too far from 1 as described in lXu et al.l (|2012l ). 
In addition, we apply a 10% Gaussian prior on 1 + e to 



Measuring Da and H using BAO 11 



limit noise from dragging e to unrealistically large values. In 
fitting the mocks without the prior, only 2 have measured e 
values > 0.1, so this prior does not have a significant impact. 
To verify that this prior is not cosmologically important, we 
use the CMB+all BAO ow w a CDM M arkov Chain Monte 
Carlo results from IMehta et all (|2012l ) to estimate the al- 
lowed distribution of e. The distribution is nearly Gaussian 
with a standard deviation of ~ 0.026 which is much less than 
our 0.1 prior. 

The monopole and quadrupole correlation function 
templates (£o,t{r) and ^2,t{r)) are derived from the 2D power 
spectrum Pt(k,fi) template (Equation [28} as described in 
32~4l We set E s = 4/iT 1 Mpc. We let f3 vary in our fits as it 
affords us leverage on the amplitude of the quadrupole with 
which it is partially degenerate. We put a prior on fj centered 
at f/b ~ fi m (z)°' 55 /6 = 0.35 before reconstruction and af- 
ter reconstruction with 0.2 standard deviation in both cases. 
The choice of /3 = as the center of the prior after recon- 
struction is because we expect the Kaiser effect to be mostly 
removed. We fix E_l = 6ft _1 Mpc and Ey = 10/i _1 Mpc in 
our pre-reconstruction fits and Ex = Em = 3/i _1 Mpc in 
our post-reconstruction fits. These values are approximated 
from the fit results to the average correlation function of the 
mocks where we set Em = (l+/)Ex in the pre-reconstruction 
case due to the Kaiser effect and Em = Ex in the post- 
reconstruction case due to the expected removal of Kaiser 
squashing by reconstruction. 

We simultaneously fit the monopole and the quadrupole 
for 4 non-linear parameters log(£>o), Pi a an d e > m addition 
to the linear nuisance parameters in Ae(r). The non-linear 
parameters are handled using a simplex algorithm and the 
linear parameters are obtained using a least-squares method 
nested within this simplex. That is, for each set of non- 
linear parameters, the least-squares algorithm returns the 
corresponding best-fit linear parameters. The simplex steps 
through the non-linear parameter space until the best-fit 
values are obtained. To determine the best-fit values, we 
minimize the x 2 goodness-of-fit indicator given by 



X = (m — d) T C 1 (rn — d) 



(52) 



where m is a column vector of the model at each step in the 
simplex and d is the data. Both of these must contain both 
the monopole and quadrupole values in sequence. C is the 
covariance matrix described below. We use a fiducial fitting 
range of 50 < r < 200/i _1 Mpc which corresponds to fitting 
50 points in both the monopole and the quadrupole. This 
gives 2 x 50 — # of fit parameters = 100 — 10 = 90 degrees- 
of-freedom (dof) in the fit. Using this technique we obtain 
best-fit values of our parameters of interest, the isotropic 
dilation a and the anisotropic warping e of the BAO signal. 

In addition, we can calculate the probability distribu- 
tion p(a, e) by fitting for the other parameters at various grid 
values of these two parameters and measuring the best-fit 
X 2 . This is feasible because p(x) oc exp(— x 2 /2). The con- 
stant of proportionality corresponds to the normalization 
that makes the integral J p(x)dx = 1. Then we can calcu- 
late p(a) and p(e) as 



p{a) 
P(e) 



p(a, e)de 
p(a, e)da. 



(53) 
(54) 



We can take the widths of these distributions (a a and cr £ ) 
as measurements of the errors on a and e if a and e have 
Gaussian posteriors. In iJ5]we demonstrate using fit results 
to our mock catalogues that a and e are consistent with hav- 
ing been drawn from Gaussian distributions. The covariance 
between a and e (Cut) can also be calculated and converted 
into a correlation coefficient p at . These are defined 



C a 



p(x)(x — {x)) 2 dx 



p(a, e)(a — — {e))dade 



(55) 
(56) 
(57) 



where (x) is the mean of the distribution p(x). For our grids 
we pick the ranges 0.7 < a < 1.3 and —0.3 < e < 0.3 at 
spacings of 0.0025 and 0.005 respectively. We also apply an 
additional Gaussian prior on log(a) with a width of 0.15 to 
suppress any unphysical downturns in the x 2 distribution 
at small a. These a correspond to the acoustic peak in the 
model being pushed out to larger scales where the fitter has 
an easier time hiding the peak in th e large errorbar s. This 
procedure was scrutinized in detail in IXu et al.l l|2012l '). Note 
that our grids are only used to compute p(a) and p(e). We 
do not use them to infer the best-fit values of these param- 
eters, which are obtained through the non-linear simplex 
optimization described above. 

We obtain a smooth estima te of the covaria nce matrix 
using the method described in IXu et all (|2012h extended 
here to include the quadrupole (see below). A detailed de- 
scription of the method is given there. This method relies 
on a form of the Gaussian covariance matrix that includes 
some additional modification parameters. These allow us to 
adjust the amount of shot-noise and sample variance to best 
match the covariances calculated directly from the mock cor- 
relation functions. To derive values for these modification 
parameters, we maximize the likelihood function 



C = n(2^" 9/2 (detC)- 1/2 exp(- X ?/2). 



(58) 



Here N is the total number of mocks and q is the number 
of points to fit. We derive the parameters using the mock 
covariances between 50 < r < 200/i~ 1 Mpc (50 monopole 
points and 50 quadrupole points) and hence q = 50 x 2 = 
100. Xi = xfC~ 1 Xi where Xi is a column vector of dimen- 
sion q containing the difference between the monopole and 
quadrupole of each mock and the average of all mocks. C is 
the modified form of the Gaussian covariance matrix derived 
below and defined in Equation (|62|l . We use this method 
since the mock covariances show evidence of noise and ide- 
ally the covariance matrix should be smooth. 

The method allows us to include the redshift depen- 
dence of the galaxy number density n assuming that it has 
no angular dependence. This is anchored in the observation 
that the covariance in configuration space is just the trans- 
form of the variance in Fourier space P 2 e :(k)/V . Hence, we 
can imagine building up the inverse variance as 



'<«-/ 



dV 
PUk) 



(59) 
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where 



dV = 



R 2 (z) 



Ho y/Q m (l + z)3 + n A 



zdzdQ. 



(60) 



for a flat universe. Here, R(z) is the comoving distance to 
redshift z. We can then redefine the variance as ^ 2 e ,(k) = 
[I 2 (k)]-\ 

Swapping this new expression for the variance into 
Equation (|44|l gives 

CiiMn),^^)) = 2(2*+l)(jtf + l) 

fc 3 dlogO) 



/ 



2ti 



(61) 



We can then insert the modification parameters cq, ci, C2 
and C3 such that 

CtfteM.&'fo)) = 2(2*+l)(2*' + l) 

fc 3 dlog(fc) 



2tt 2 



Ji(kri)J e ,(k rj 



• < P«'(fc;co,ci,c 2 ) 



+ c 3 . 



(62) 



Here, where we have made the substitution 



P(k,p) + - 
n 



CoPdw(k, jit) + -^-r 



+ 



C2 

n(z) 



(l + ^ 2 ) 2 F(fc,/i,S s ) 



(63) 



in Equation ()34p . The Co term adjusts the magnitude of the 
sample variance. The ci term acts like a non-linear shot- 
noise component and the C2 term adjusts the magnitude of 
the standard Poisson shot-noise contribution. The C3 term 
can be associated with the integral constraint, which appears 
as an additive offset in the correlation function. 

In calculating the covariance matrix, we set j3 = f/b 
before reconstruction and j3 = after reconstruction, again 
due to the expected removal of large-scale redshift-space dis- 
tortions. We test several cases where we vary ft from these 
fiducial values and find that changing /3 affects the relative 
amplitudes of the Cjj(£o(ri),£oOj)) and C i3 ■(^(n), &(rj)) 
terms. This can cause slight changes in the resulting a a and 
cr e at the 0.1% level, which is not significant at our current 
levels of statistical precision. The E s streaming scale for the 
FoG is fixed at 4h~ Mpc. We find very little difference in 
the resulting modification parameters in cases where we al- 
low E s to vary. We fix Ex and E11 to their fiducial model 
values (recapped below) in our covariance matrix calcula- 
tions. The modification parameters we obtain are Co = 1.06, 
ci = 0.11, C2 = 1.49, C3 = 5.18 x 10 -8 before reconstruction, 
and c = 1.12, ci = 0.05, c 2 = 1.58, c 3 = 8.82 x 10" 8 after 
reconstruction. With these modification parameters in hand, 
we can construct a smooth approximation to the mock co- 
variances from the binned Gaussian covariance matrix using 
Equation ((621) . 

Our fiducial model parameters are summarized as fol- 
lows: we define our fiducial model before reconstruction to 
have E x = 6/i -1 Mpc, E y = lO/i^Mpc, E s = 4/i _1 Mpc and 
a prior on P centered on /? = 0.35. After reconstruction, we 
set Ex = Em = 3/i -1 Mpc and center our j3 prior on 0. Our 
fiducial fitting range is 50 < r < 200/i _1 Mpc. 

In Sj5] and SjfJ we present our measured values of a and 



e for the mocks ( jj4.1|) and actual survey data ( i|4.2p using 
the fitting models defined in Equations (|49|) & (|50|) and the 
modified Gaussian covariance matrix described above. 



4 DATASETS 
4.1 Simulations 

We use the Large Suite of Dark Ma tter Simulations (Las- 
Damas; iMcBride et al.l 12012. in pred ) to calibrate and test 
our reconstruction parameters, fitting template and covari- 
ance matrix. The LasDamas collaboration Q has provided 
publicly available mock galaxy catalogues based on these 
simulations for the Sloan Digital Sky Survey (SDSS) data 
release 7 (DR7) luminous red galaxy (LRG) sample. In par- 
ticular, we use the "gamma" release mock catalogues re- 
ferred to as lrgFull. 

The LasDamas simulations were run assuming a flat 
ACDM cosmology with fl b = 0.04, fi m = 0.25, h = 0.7, 
n 3 — 1.0 and as = 0.8. Although various box sizes were 
implemented, the 40 simulations used to construct the LRG 
mocks were 2.4/i _1 Gpc on a side with 1280 3 particles in 
each. The initial particle positions were set using second- 
order Lagrangian perturbation theory at z — 49. To con- 
struct mock galaxy catalogues from the simulations, the 
dark matter halos were populated according to halo occu- 
pation parameters tuned to match the observed clustering 
of the DR7 LRGs. In addition, the mock catalogues include 
observational effects such as redshift-space distortions and 
mimic the angular selection function of the LRG sample. 
The redshift range covered by the mocks is 0.16 < 2 < 0.44. 
We note that this is slightly different to the flux-limited 
LRG sample described in the following section that will be 
employed in this study. When fitting the data, we account 
for this by extrapolating our covariance matrix using the 
observed n(z) and the formalism described in the previous 
section. In addition, we slightly downsample the n(z) dis- 
tribution in the mocks to better match the data. Our re- 
gion of interest, the SDSS Northern Galactic Cap, covers 
~ 7200 deg 2 on the sky. The resulting geometry allows 4 
mocks to be constructed from each simulation and hence we 
have a total of 160 mocks for our analysis. 



4.2 SDSS DR7 

The observationa l data set used in this study is the SDSS 
(|York et alj200Ch DR7 (|Abazaiian et al.ll2009r ) LRG sample. 
The same dataset was also used i n the monopole-only B AO 
analysis of IPadmanabhan et al.l (|2012T l. IXu et al.1 ( |2012h & 
iMehta et al.l (|2012T ). 

The SDSS has taken photometric observations of ~ 
10, 000 deg 2 on the sky and obtained spectroscopic fol- 
lowup of nearly a million of these detected obj ects. It uses 
a dedicated 2.5m telescope (|Gunn et al.l l2006h at Apache 
Point O bservatory which h as a specially designed wide field 
camera (IGunn et al. 1 119981 ). Photometric observations were 
taken in the ugriz bands (|Fukugita et al.|[l996l ; ISmith et al.l 
120021) b y drift scanning the sky under favourable con- 
ditions (|Hogg et al.1 l200lf ). These images were then fed 



http: / /lss.phy. vanderbilt.edu/lasdamas 
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through an automated pipeline that performed the neces- 
sary astrometric and photometric calibrations. The pipeline 
also detected and measured the photometric properties of 



also detected and mea sured the p hotometric proper ties or 
the observed objects JPier et all 120031; llvezic et all 12004 
iTucker et al.| |2006l; [Padmanabhan et al l 120081'). Selec t sub- 
samples |Strauss et al.l |2002| ; lEisenstein et al.l l200ll ) were 
then designated for spectroscopic followup using a 640 fiber 
spectrograph. 

The DR7 LRG sample is part of the last data release of 
SDSS-II, the second phase of SDSS which was completed in 
2009. The LRG sample was select ed according to the pre- 
scription in lEisenstein et ah! (|200ll ). This selection was opti- 
mized to identify the most luminous (and hence most mas- 
sive and highly biased) galaxies which can be observed out 
to high redshifts. Since the volume encompassed by equal 
angles on the sky increases with redshift, we can probe the 
large volumes necessary for cosmological studies using these 
luminous galaxies. The LRGs tend to be old systems with 
uniform spectral energy distributions that exhibit a strong 
4000A break. This gives them a distinct colour-flux-redshift 
relation which allows them to be uniformly selected over 
a wide redshift ran ge. Our sample matches exactly that of 
iKazin et all ^OlOal ) and we refer the interested reader there 
for details of its construction. We use the flux-limited LRG 
sample in the SDSS Northern Galactic Cap only. This sam- 
ple spans a redshift range of 0.16 < z < 0.47 and has a 
number density of ~ 10 -4 /i 3 Mpc -3 . 



5 MOCK CATALOGUE RESULTS 

In this section we present the results of fits to the LasDamas 
mock correlation functions before and after reconstruction. 
These were computed and fit by taking the LasDamas cos- 
mology as the fiducial cosmology. 

We plot the average monopole and quadrupole of the 
160 mocks before reconstruction in Figure [4] The fiducial 



templates (solid lines in Figures 1(b) & [2]) are shown as 



the grey dotted lines. The monopole and the quadrupole at 
large scales look similar to the fiducial templates. However, 
the small scales (r < 50/i~ 1 Mpc) in the quadrupole show 
substantially different structure indicating that our model 
does not fit the data well at these scales. This motivates our 
choice for the fiducial fitting range : 50 < r < 200/i~ 1 Mpc. 
The best-fit model to the monopole and quadrupole are over- 
plotted as the red lines where the solid line corresponds to 
using the fiducial Ao t 2(r) and the dashed line corresponds to 
using Ao,2(r) = (i.e. no broadband marginalization). The 
monopole fits look very similar; however the fiducial Ao : 2(r) 
does much better in the quadrupole near the acoustic scale 
with x 2 decreasing by ~ 33 relative to the Ao,2{r) = fit. We 
allow Ex and E11 to vary in these fits and obtain best-fit val- 
ues of 6.3h~ Mpc and 10.4/i _1 Mpc, motivating our choices 
in the fiducial model. We keep E s fixed at 4h~ Mpc. 

A comparison of the monopole and quadrupole be- 
fore and after recon s truct ion is shown in Figure [5] As in 
iPadmanabhan et al.l (|2012h . we see the acoustic peak in the 
monopole appears less smeared after reconstruction which 
indicates that our reconstruction technique was effective 
at partially undoing non-linear evolution. This is also re- 
flected in the fact that after reconstruction (where we as- 
sume the smearing is isotropic), a fit to the average of 



the mocks gave a much smaller BAO smoothing scale of 
E n i = 2.9fr _1 Mpc as opposed to the pre-reconstruction val- 
ues of E ± = 6.3/i _1 Mpc and Ey = 10.4/i _1 Mpc. In addition, 
we see that the quadrupole is nearly zero on large scales 
after reconstruction, which implies that our reconstruction 
technique was also effective at partially removing the Kaiser 
effect. The fact that the quadrupole is positive and not ex- 
actly is likely due to some slight anisotropy introduced by 
reconstruction. This is discussed in more detail below and 
shown not to significantly affect our measurements of a and 
c. 

By averaging the mocks we have effectively increased 
the survey volume by a factor of 160 and hence the variance 
should decrease by an equal amount. This means that the 
average of the mocks should have substantially less noise. 
In addition, we know that we are computing the correlation 
functions and fitting using the correct (LasDamas) cosmol- 
ogy. Hence, we expect that the a and e values measured from 
the average of the mocks should be 1 and respectively if our 
models are unbiased (i.e. there should be no shift in the lo- 
cation of the peak relative to the model and there should be 
no anisotropy). We find that fitting the pre-reconstruction 
mock average gives a = 1.005 and e = 0.003 while post- 
reconstruction we measure a — 1.002 and e — 0.002. The 
slight offset in a from 1 is not too concerning as we expect 
non-linear structure growth to shift the peak b y < 0.5% 
(|Padmanabhan fc Whitell2009] ; iMehta et al.ll201lh . The fact 
that a moves closer to 1 after reconstruction is encourag- 
ing as reconstruction is supposed to remove some effects of 
non-linear structure growth. The small bias in e is not sig- 
nificant at our current levels of statistical precision and is 
likely the result of small mismatches between the broadband 
model and the data. This will be discussed in more detail 
shortly. We again emphasize that in these fits we have al- 
lowed E|| and Ex to vary. When we fit the individual mocks, 
the signal-to-noise of the data is not sufficient for constrain- 
ing any of these parameters and hence we fix them in the 
fiducial model to the values obtained in the averaged mock 
fits. 

We measure a and e for each mock using the fitting pro- 
cedure and fiducial model outlined in £13.31 We also estimate 
the errors a a , a t and the correlation coefficient p at for each 
mock using Equations (|55|) & (|5T[) . Before reconstruction, 
we measure a mean (a) = 1.003 ± 0.003 with an rms scatter 
between the mocks of 0.034 and a median a — 1.008 with 
16th/84th percentiles of the mocks corresponding to 1q 1)36 
(these will henceforth be denoted the quantiles). For e we 
measure a mean (e) = 0.001 ± 0.003 with an rms scatter 
between the mocks of 0.037 and a median e = 0.004 with 
quantiles to 037. After reconstruction, we measure a mean 
(a) = 1.002 ± 0.002 with an rms scatter between the mocks 
of 0.024 and a median a = 1.002 with quantiles to.oll- For e 
we measure a mean (e) = 0.002 ± 0.003 with an rms scatter 
between the mocks of 0.032 and a median e = 0.007 with 
quantiles to'ca?- These values were calculated after rejecting 
the mocks where the acoustic signal is too weak to obtain an 
accurate centroiding of the BAO peak. This corresponds to 
making a cut in a a at 0.07 and disc arding the mock s that lie 
above this cut as demonstrated in IXu et al.l (|2012l ). Before 
reconstruction, 9 mocks out of 160 lie above this cut and 
after reconstruction, lie above this cut. 

The median e in the pre- and post-reconstruction cases 
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Figure 4. Average monopole (left) and quadmpolc (right) of the 160 mock catalogues before reconstruction. The monopole and the 
quadrupolc at large scales are similar to the fiducial templates (grey dotted lines, identical to the solid lines plotted in Figurcs |l(b)| fc 
The quadrupole on small scales (r < 50h Mpc), however, shows substantially different structure to the fiducial template. The fit to 
the average of the monopole and quadrupole from the mocks is overplotted in red. The solid line corresponds to a fit using the fiducial 
•Ao,2(r) (Equation H51I and the dashed line corresponds to a fit using Ao,2{r) = 0. We allow Ex and Sii to vary in these fits and obtain 
best-fit values of 6.3h~ Mpc and 10Ah~ Mpc respectively using the fiducial Ao,2(0- In the monopole case, the fit using the fiducial 



Ao,2(t) is very similar to the Ao^ir) = 
decreased by ~ 33 going from Aq,2(0 



fit. In the quadrupole, the Ao,2(r) 
= to the fiducial Ao,2(f)- 



fit is much worse around the acoustic scale. Overall, \ 
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Figure 5. The average monopole (left) and quadrupole (right) of the 160 mocks before (gray crosses) and after (black crosses) recon- 
struction. One can see that after reconstruction, the acoustic peak in the monopole has sharpened up, indicating that reconstruction is 
effective at removing the degradation of the BAO caused by non-linear structure growth. In the quadrupole, the power at large-scales 
goes close to which implies that reconstruction was effective at removing the Kaiser effect. It is not exactly zero due to some small 
anisotropy introduced by the reconstruction technique itself (sec Figure |BJ. We note that the quadrupole is multiplied by r 2 in this figure 
and hence the magnitude of this anisotropy is exaggerated. 



are different from and from the mean at > 1-2 times the er- 
ror on the mean. In addition, the post-reconstruction quan- 
tiles are asymmetric, implying that the posterior e distri- 
bution deviates from Gaussian. These appear to be in part 
due to the intrinsic noise in the data and in part due to a 
slight mismatch between the model and the data. To reduce 
noise, we effectively increase the spatial volume of the data 
by combining our 160 mocks into groups of 2, 4 and 8, and 
re-perform our fits. In general, we see a better agreement be- 
tween the mean and median e. The quantiles remain mildly 
asymmetric in some cases but overall we see improved agree- 
ment. The rms scatter decreases by roughly the expected 



amount (~ \J~rn, where m = 2, 4 or 8) if we consider e to be 
Gaussian. These results are summarized in Table [T] 

We see that there is a persistent bias in e towards non- 
zero values that is currently below our detection thresh- 
old. This bias is < la significant before reconstruction and 
only at the 1-1. 5a level after reconstruction. To further test 
this, we split the mocks into 2 groups of 80 which reduces 
the noise in the data. After re-performing the fits, we find 
(e) ~ 0.002 both before and after reconstruction. These val- 
ues agree with the fit results to the average of the 160 mocks 
described above. This suggests that the persistent bias in e is 
not due to noise but is rather a result of some mismatch be- 
tween the model and the data. In our fits we fix Ex, Ey and 
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Figure 6. Average quadrupole from the mocks in real space be- 
fore (grey) and after (black) reconstruction. The quadrupole be- 
fore reconstruction is very close to as we would expect in real 
space due to the lack of redshift-space distortions. Our e measure- 
ments are unbiased in this case which suggests that the small bi- 
ases we see in rcdshift space are due to slight mismatches between 
our redshift-space distortion models and the actual broadband in 
the data. After reconstruction, the quadrupole at large scales ac- 
quires some additional power likely due to the survey geometry 
and sample number density fluctuations as a function of redshift. 
Our post-reconstruction real space e measurements remain unbi- 
ased which suggests that this anisotropy can be accounted for by 
our A2(r) nuisance terms. 




Figure 7. /3 versus e measured from fits to the 160 mocks after 
reconstruction. One can see that these two parameters are un- 
correlated. We see similar results before reconstruction as well. 



E s , and use the nuisance terms in A2(r) to account for any 
other mismatch in the broadband shape between the model 
and the data. However Ex, E|| and E s are partially degen- 
erate with e, so if they are fixed at non-optimal values that 



cannot be fully compensated by A2(r) (see Figure 12(b) I, 
the fitter can resort to adjusting e. We stress however, that 
such small biases in our redshift-space measurements of e are 
below the current detection limit in a single DR7 realization 
as indicated by the rms of the mocks. 

We can gain additional insights by analyzing the real- 
space mocks which do not have redshift-space distortions 
and therefore do not require E s or anisotropic E n i in the 
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Figure 8. a versus e (top) and Da{z) versus H(z) (bottom) for 
the mocks after reconstruction. One can see that a and e are not 
highly correlated. From the points plotted we measure a correla- 
tion coefficient of 0.27 and analogously in the prc-reconstruction 
case we measure 0.20. These are in excellent agreement with 
Fisher matrix predictions (p at ~ 0.21). We see a stronger cor- 
relation between D \ and H which we obtained by combining a 
and t as in Equations J5} & ©. We expect Pd a h ~ 0.4 and we 
find correlation coefficients of 0.23 and 0.50 between our measured 
D \ and H values before and after reconstruction respectively. 



Table 1. e statistics for various mock combinations. The first 
column indicates the number of mocks we have combined (m). 
The second column quotes the mean e we measure with the stan- 
dard error on the mean. The third column shows the rms of the 
mocks. The fourth column quotes the median e and the fifth col- 
umn quotes the quantilcs. 



Qtls 



Redshift Space without Reconstruction 



1 


0.001 ± 0.003 


0.037 


0.004 


+0.032 
-0.037 


2 


0.001 ± 0.003 


0.029 


0.006 


+0.023 
-0.030 


4 


0.001 ± 0.003 


0.019 


-0.001 


+ 0.017 
-0.010 


8 


0.002 ± 0.003 


0.012 


0.002 


+0.009 
-0.013 



Redshift Space with Reconstruction 



0.002 ± 0.003 
0.003 ± 0.002 
0.003 ± 0.002 
0.003 ± 0.002 



0.032 
0.018 
0.012 
0.008 



0.007 
0.006 
0.003 
0.006 



+0.023 
-0.037 
+ 0.012 
-0.020 
+0.012 
-0.011 
+0.004 
-0.007 
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Figure 9. Post-reconstruction a a versus cr e for the mocks. We see 
that the errors on a and the errors on e are directly correlated with 
each other. This indicates that mocks with poorer measurements 
of the acoustic scale a also have poorer measurements of the BAO 
anisotropy e. A similar correlation exists in the pre-reconstruction 
case. The median ratios of o" e /(l + e)-to-cr a /a are ~ 1.3 and 
~ 1.4 before and after reconstruction respectively. Fisher matrix 
arguments predict a ratio of ~ 1.2. 



model. We find that these give non-biased measures of e 
in both the pre- and post-reconstruction cases. In the pre- 
reconstruction case we measure the mean e to be (e) = 
0.003 ± 0.003 with an rms between the mocks of 0.037. The 
median is e = 0.002 with quantiles Iq'oIo- After reconstruc- 
tion we measure (e) = 0.001T0.002 with a mock rms of 0.027 
and e = —0.002 with quantiles 023- One can see that the 
mean and median e are consistent with each other and with 
0. Fitting the average of the 160 mocks gives e = 0.001 and 
0.000 before and after reconstruction respectively; again im- 
plying a largely unbiased measurement of e in real space. 

An interesting artifact we do find is that reconstruction 
appears to introduce some broadband anisotropy as shown 
in Figure [6] Here we have plotted the mean of the real-space 
quadrupoles before (grey) and after (black) reconstruction. 
We see that the quadrupole is nearly before reconstruc- 
tion as expected since there should not be any anisotropies 
in real space. However, after reconstruction, the quadrupole 
acquires some additional large-scale power. The reconstruc- 
tion displacement vectors may take on a subtle anisotropy 
due to the variation of number density with redshift or the 
survey geometry (i.e. if it is wider than it is deep). For- 
tunately, this broadband anisotropy introduced by recon- 
struction is fairly smooth and can be removed by the A2(r) 
nuisance parameters as evidenced by our unbiased measures 
of e in post-reconstruction real space. 

To build more intuition for the parameters we fit in 
redshift space and to demonstrate their inter-dependencies, 
we show various scatter plots of these quantities in Figures 
□ H El HQ] & E] Figure [7] shows the values of /3 and e we 
obtain from our fits to the 160 mocks after reconstruction. 
Our pre-reconstruction results are similar. One can see that 
these two parameters are not correlated with each other. 

The top panel of Figure[8]shows the e versus a values we 
measure from the mocks after reconstruction. Again we see 



that these two parameters are not highly correlated. The 
correlation coe fficient between Da and H is predicted to 
be pd a h ~ 0.4 l)Seo fc Eisensteinll2007l ). This subsequently 
predicts a oh / 'H-to-an A / Da ratio ~ 2 (i.e. the fractional 
error of the Hubble parameter is twice that of the angu- 
lar diameter distance). Using these values, a Fisher matrix 
argument shows that we should expect p at ~ 0.21 (see Ap- 
pendix The correlation coefficient between the a and e 
values plotted in Figure [5] is 0.20 and the corresponding pre- 
reconstruction value is 0.27. Both are in excellent agreement 
with the Fisher matrix prediction. The bottom panel of Fig- 
ure [S] shows our a and e measurements translated into mea- 
surements of Da and H using Equations ([5]) & ((6j). In the 
plotted post-reconstruction case, the correlation coefficient 
between Da and H is ~ 0.50 and in the pre-reconstruction 
case it is ~ 0.23, which are not too different from our as- 
sumed Pd a h = 0.4. 

Figure [9] shows a e versus a a for the mocks after recon- 
struction. We see that the errors on a and e are correlated 
which implies that mocks with poorer measurements of the 
acoustic scale (i.e. larger <j a values) also have poorer mea- 
surements of the BAO anisotropy (i.e. larger a t values). We 
see a similar correlation in the pre-reconstruction results. 
Taking the ratio of c e /(l + e)-to-<r Q /a, we find a median 
~ 1.3 before reconstruction and ~ 1.4 after reconstruction. 
Fisher matrix arguments predict a ratio of ~ 1.2 (Equation 
IA13|I . which is similar to what we see. 

Figure [TO] shows the values of e we measure versus a e 
before (left) and after (right) reconstruction. Reconstruc- 
tion clearly decreases the scatter in e which is again high- 
lighted in Figure [TT] showing the values of e before and af- 
ter reconstruction. While the two are correlated, the post- 
reconstruction values have smaller scatter as evidenced by 
the locus of points having a slope shallower than 1:1. From 
Figure [10] we also see that reconstruction decreases cr £ , our 
estimated error on e. However, these a e values are still large 
compared to a a (see Figure 0. 

Next we test the robustness of our fitting model to 
changes in various model parameters. A full list of these 
results are found in Table [2] for changes of Ex, £11, Ss, 
fitting range and form of Ai(r) both before and after re- 
construction. In the table, poly2 corresponds to an Ai{r) — 
ai/r 2 + a-ijr and polypi corresponds to an A%(r) = ai/r 2 + 
a.1 1 r + (13 + a^r. We see that the scatter in e between the 
mocks can show ~ 10% variations quite often. This again 
indicates the noisiness of our e measurements. 

Figure [T2l shows scatter plots in a and e for a few sample 
cases. Here, Aa and Ae are the differences between the a 
and e values measured using the slightly altered model and 
the fiducial model. We expect the average Aa and Ae to 
be within the errors if our measurements of a and e are 
consistent between the two models. We see that in all cases, 
Aa = and Ae = fall within the scatter predicted by 
the quantiles on a mock-by-mock basis. However, the errors 
on the average Aa and Ae are on the order of yl60 times 
smaller than the scatter implied by the mocks. This indicates 
that in a few cases, we detect a significant shift in the average 
value of a and e measured. 

In particular, this occurs in the pre-reconstruction cases 
where we have changed the fitting range. The a and e scat- 
ter plots for the 30 < r < 200/i -1 Mpc fitting range versus 



the fiducial fitting range cases are shown in Figure 12(a) 
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Figure 10. e versus a t before (left) and after (right) reconstruction for the 160 mocks. One can see that reconstruction decreses the 
scatter in the measured t values (this is further highlighted in Figure [TTJ. While reconstruction does decrease the average error on e, we 
see that the errors we measure are still fairly large compared to the errors on a (see Figure [9}. 
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Figure 11. e before reconstruction versus e after reconstruction 
fit from the 160 mocks. The slope of a linear fit to these points 
is less than 1 implying that the post-reconstruction e values have 
smaller scatter. 



The plots are shown in pairs with a on the left and e on 
the right. The first plot in each pair corresponds to the 
pre-reconstruction case and the second plot to the post- 
reconstruction case. We see that on average, the larger fit- 
ting range gives slightly smaller values of a and e. If we 
begin fitting at r = 30/i -1 Mpc where the errorbars are 
smaller, the fitter forces the model to match the data at 
these small scales where we know the templates (especially 
the quadrupole) are not faithful representations of the data. 
The Ao^ir) marginalization terms compensate for this at 
the expense of accurately fitting the BAO scale. After re- 
construction, Aa and Ae are both which suggests that 
the model is better matched to the data in this case. 

We also see average Ae values that are significantly 
different from zero in the E s = 0/i _1 Mpc case both be- 
fore and after reconstruction. This is illustrated in Figure 
= 0/i -1 Mpc implies that we exclude FoG from the 



The likely culprit here is again the mismatch between the 
data and the model at small r especially in the quadrupole. 
In addition, if we compare the dotted line in Figure |2(b)| 
(the E s = 0/i _1 Mpc case) and the average quadrupole in 
Figure U we see that the quadrupole BAO feature in this 
model is a poorer fit to the data overall. The mocks show 
more of a crest-trough-crest structure near the BAO scale 
as in the fiducial parameter template (solid lines in Figure 
[2]| whereas the trough in the E s = 0/i -1 Mpc case is much 
weaker. This is further affirmed by the fact that the dot- 



ted line in Figure 2(a) (isotropic E n i) has a similar looking 



turn 



model, which is unrealistic as it is implemented in the mocks. 



BAO feature and shows a similar discrepancy in e relative 
to the fiducial model before reconstruction. The fitter can 
partially compensate for these differences through adjusting 
the value of e which also gives rise to crests and troughs 
near the BAO scale, although with different structure than 
those introduced through Ex, Ey and E s . Hence one must 
pick a quadrupole model that has a BAO feature fairly well 
matched to the data to avoid biasing the e values measured. 

Despite these offsets in the median a and e for differ- 
ent fitting models, we note that at the statistical precision 
of current datasets, we would not be able to detect any of 
these changes. For the DR7 mocks, the average a a is ~ 0.04 
before reconstruction and ~ 0.03 after reconstruction. The 
average a t are even larger at ~ 0.05 before reconstruction 
and ~ 0.04 after reconstruction. Hence, assuming that a a 
and cr e characterize the error on a and e, a 0.003 shift will 
fall entirely within the expected errors. Therefore, our fitting 
model is reasonably robust against small changes to model 
parameters and our measured a and e values are largely 
unbiased. 

We perform similar exercises for various different re- 
construction parameters such as the bias and (3 values we 
input to the algorithm. The fiducial reconstruction param- 
eters we use are b — 2.2 and j3 — 0.3. We also generate the 
constrained Gaussian realizations in a manner equivalent 
to Wiener filtering the data in the observed region. This is 
akin to setting the unconstr ained Gaussian realization am- 
plitudes (8u in Equation 3 of lPadmanabhan et al.l (|2012h 'l to 
0, the expected mean amplitude from averaging over many 
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Table 2. Fitting results from the mocks for various models. The model is given in column 1. The median a is given in column 2 with 
the 16th/84th percentiles from the mocks given in column 3 (these are denoted as the quantiles in the text, hence the label Qtls in the 
tabic). The median e is given in column 6 with corresponding quantiles in column 7. The median difference in a on a mock- by-mock 
basis between the model listed in column 1 and the fiducial model is given in column 4 with corresponding quantiles in column 5. The 
analogues for e are given in columns 8 and 9. The mean x 2 /dof is given in column 10. 



Model 




a 


Qtls 


Aa 


Qtls 


€ 


Qtls 


Ae 


Qtls 


(x 2 )/dof 


Redshift Space without Reconstruction 


Fiducial [/] 




1.008 


+0.030 
-0.036 






0.004 


+0.032 
-0.037 






91.92/90 


Fit w/ (E±,S||) -> (8, Sjft" 1 


Mpc. 


1.008 


+0.029 
-0.039 


0.001 


+0.003 
-0.003 


0.001 


+0.032 
-0.037 


-0.002 


+0.003 
-0.003 


91.97/90 


Fit w/ S s -> O/i^Mpc. 




1.005 


+0.030 
-0.037 


-0.002 


+0.003 
-0.004 


0.000 


+0.031 
-0.033 


-0.003 


+0.005 
-0.004 


91.76/90 


Fit w/ A2(r) = poly2. 




1.006 


+0.031 
-0.036 


0.000 


+0.002 
-0.001 


0.005 


+0.032 
-0.036 


0.001 


+0.005 
-0.005 


93.08/91 


Fit w/ A2(r) = polyA. 




1.006 


+0.031 
-0.035 


0.000 


+0.001 
-0.001 


0.002 


+0.037 
-0.033 


-0.000 


+0.007 
-0.007 


90.98/89 


Fit w/ 30 < r < 200ft" 1 Mpc 


range. 


1.003 


+0.032 
-0.037 


-0.003 


+0.004 
-0.005 


0.000 


+0.032 
-0.035 


-0.003 


+0.005 
-0.005 


105.90/104 


Fit w/ 70 < r < 200ft" 1 Mpc 


range. 


1.007 


+0.030 
-0.035 


0.000 


+0.002 
-0.002 


0.004 


+0.032 
-0.040 


0.000 


+0.003 
-0.003 


79.33/76 


Fit w/ 50 < r < 150ft _1 Mpc 


range. 


1.005 


+0.031 
-0.040 


-0.000 


+0.004 
-0.006 


0.004 


+0.036 
-0.043 


-0.001 


+0.008 
-0.007 


54.34/58 


Redshift Space with Reconstruction 


Fiducial [/] 




1.002 


+0.023 
-0.022 






0.007 


+0.023 
-0.037 






92.68/90 


Fit w/ (S ± ,E||) -> (2,4)h~ 1 


Mpc. 


1.002 


+0.023 
-0.021 


-0.000 


+0.001 
-0.001 


0.008 


+0.022 
-0.037 


0.000 


+0.001 
-0.001 


92.73/90 


Fit w/ S s -> Oft^Mpc. 




1.001 


+0.024 
-0.020 


-0.001 


+0.005 
-0.003 


0.004 


+0.020 
-0.030 


-0.003 


+0.007 
-0.005 


92.27/90 


Fit w/ A2(r) = poly2. 




1.002 


+0.024 
-0.022 


0.000 


+0.001 
-0.000 


0.008 


+0.023 
-0.037 


0.001 


+0.002 
-0.002 


94.23/91 


Fit w/ A2(r) = polyA. 




1.002 


+0.023 
-0.022 


0.000 


+0.001 
-0.001 


0.005 


+0.025 
-0.037 


-0.001 


+0.004 
-0.004 


91.68/89 


Fit w/ 30 < r < 200ft" 1 Mpc 


range. 


1.003 


+0.022 
-0.023 


0.000 


+0.002 
-0.001 


0.006 


+0.022 
-0.038 


0.000 


+0.002 
-0.002 


106.12/104 


Fit w/ 70 < r < 200ft" 1 Mpc 


range. 


1.002 


+0.023 
-0.021 


0.000 


+0.001 
-0.001 


0.006 


+0.022 
-0.036 


-0.000 


+0.002 
-0.002 


79.99/76 


Fit w/ 50 < r < lSOft^Mpc 


range. 


1.002 


+0.024 
-0.023 


-0.001 


+0.004 
-0.004 


0.008 


+0.022 
-0.037 


-0.001 


+0.005 
-0.005 


54.65/58 


Recon. w/ -> 0.24. 




1.002 


+0.023 
-0.022 


-0.000 


+0.001 
-0.001 


0.005 


+0.023 
-0.033 


0.000 


+0.002 
-0.003 


92.49/90 


Recon. w/ -> 0.36. 




1.002 


+0.022 
-0.020 


0.000 


+0.001 
-0.002 


0.005 


+0.024 
-0.038 


-0.000 


+0.002 
-0.002 


92.89/90 


Recon. w/ b — > 1.8. 




1.001 


+0.022 
-0.021 


-0.000 


+0.006 
-0.005 


0.006 


+0.025 
-0.041 


-0.000 


+0.006 
-0.006 


92.61/90 


Recon. w/ b -+ 2.6. 




1.003 


+0.023 
-0.023 


0.001 


+0.004 
-0.004 


0.004 


+0.025 
-0.036 


-0.001 


+0.006 
-0.005 


92.65/90 


Recon. w/ Wiener Filter. 




1.004 


+0.020 
-0.022 


-0.000 


+0.004 
-0.003 


0.005 


+0.024 
-0.035 


-0.000 


+0.004 
-0.004 


92.69/90 


Recon. on f! m = 0.4 case. 1 




0.832 


+0.021 
-0.019 


-0.171 


+0.010 
-0.010 


0.020 


+0.028 
-0.037 


0.017 


+0.014 
-0.014 


92.61/90 



a = 1 and e = in the LasDamas cosmology correspond to a = 0.832 and e = 0.013 in this f2 m = 0.4 cosmology according to 
Equations J5]&[6jl. 



constrained realizations. In our fiducial procedure, the Su 
are matched to the observed density field. The fitting re- 
sults of these cases are also listed in Table We see that 
in all cases we recover similar values of a and e. This indi- 
cates that our fitting model effectively marginalizes away the 
broadband signal introduced by reconstruction when incor- 
rect values of the fiducial parameters are used. The model is 
robust against small changes in reconstruction parameters 
and returns reliable measurements of a and e. 

Finally we calculate and perform our fits using a fiducial 
cosmology that is significantly different to the LasDamas 
cosmology. This forces a stronger anisotropic BAO signal 
to appear in the quadrupole. We pick a cosmology with 
fim = 0.4 that preserves the matter-to-baryon ratio of Las- 
Damas. We also fix Q. m h 2 which implies h = 0.553. We 



convert the measured a and e values to Da{z) and H(z) 
using Equations ([5} & © and compare these to the values 
measured using the fiducial cosmology. This is illustrated in 
Figure [131 The equations listed above only allow us to infer 
DA(z)/r s and H(z)r s . We have assumed r s = 159.71 Mpc, 
which is the sound horizon in the LasDamas cosmology, to 
obtain the Da{z) and H(z) values plotted in the figure. In 
the LasDamas cosmology (which is the true cosmology in our 
mocks), D A {z) = 1032 Mpc and H(z) = 81.8 km/s/Mpc at 
z = 0.35. Taking the ratio AL)2/Da(z) and AH/H(z) im- 
plies that on average our measurements of Da (z) and H(z) 
using the Sl m = 0.4 cosmology and the true LasDamas cos- 
mology differ by ~ 0.7% and ~ 0.6% respectively. Dividing 
the scatter indicated by the quantiles by V160 suggests that 
these average offsets are significant, although again, in a sin- 
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(a) Comparison of results obtained using a larger fitting range (30 < r < 200/i 
average Aa that is ~ —0.003 before reconstruction. The error on the mean is a factor of ~ %/160 smaller than the scatter indicated by 
the quantiles, which makes this 0.003 shift statistically significant. This is caused by the templates being poor matches to the data at 
low r. While the Aq 2(1") terms attempt to compensate for this, accurate fitting of the BAO scale is compromised. The significant Ae 
in this case is rooted in the same cause. The fact that Aa and Ae are both after reconstruction suggests that the post-reconstruction 
model is better matched to the data. 
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(b) Comparison of results obtained using S s = O/i Mpc versus the fiducial value of S s = Ah~ Mpc. Here, we see that the average Ae is 
different from when considering the error on the mean which is a factor of ~ V160 smaller than the scatter indicated by the quantiles. 
This is caused by the quadrupole model being a less optimal match to the data. We note that these small shifts in a and e are certainly 
not detectable in each mock, which have much larger errors on a and e than ~ 0.003. 



Figure 12. The robustness of our fitting model as demonstrated by a and e scatter plots for some sample cases. Results for other 
changes to the fiducial model are given in Table f2] The plotted a and e values were obtained through fitting the mock correlation 
functions before and after reconstruction. The plots are presented in pairs: a on the left and e on the right. The first plot in each pair 
shows the pre-reconstruction results and the second plot shows the post-reconstruction results. The black crosses indicate the medians 
and quantiles of the mock measurements. If we obtain consistent measurements of a and e with a model that has parameters slightly 
different to the fiducial model, then we should see Aa and Ae values that are ~ 0. We see that this is true at the level of our current 
statistical precision. Therefore, our fitting model is reasonably robust against small changes to the fiducial model parameters. 



gle mock, we would not be able to detect these offsets. The 
fiducial Da(z) and H(z) are calculated assuming a median 
redshift of z = 0.35; however, if the true median redshift 
were slightly different, such discrepancies would not be un- 
expected. In addition, our models for a and e are based on 
Taylor expansions around 1 and respectively. When the 
fitting model is constructed using a fiducial cosmology that 
is extremely wrong, the a and e values we measure will de- 
viate significantly from 1 and as shown in the last row of 
Table [21 Our first-order Taylor assumption may be break- 
ing down at this point, further affecting our measurements. 
In this case, one can iteratively change the fiducial cosmol- 
ogy and re-fit for a and e until values closer to 1 and are 
obtained. 

We verify our assumption that the second moments of 
p(a) and p(e) are good indicators of the errors on a and e in 
Figure 1141 The top panel of this figure shows a normalized 



histogram of (a — {a))/a a after reconstruction and the bot- 
tom panel shows the analogue for e. The best-fit Gaussian is 
overplotted. We perform K-S tests on these distributions and 
list the p- values in the plots. These give the probability that 
the plotted distribution is drawn from a Gaussian. One can 
see that the post-reconstruction p-values are 0.50 and 0.20 
for a and e respectively. For comparison, the corresponding 
pre-reconstruction p-values are 0.40 and 0.09. These values 
indicate that there are finite probabilities that a and e have 
Gaussian posteriors. Hence the standard deviations a a and 
cr e we calculate from x 2 (o;,e) characterize the errors on a 
and e reasonably well. 
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Figure 13. D^(z) and H(z) scatter plots obtained by plotting 
those measured using an Q m = 0.4 cosmology versus the true 
LasDamas cosmology. Da{z) is in units of Mpc and H(z) is in 
units of km/s/Mpc. These values were calculated using Equations 
© & © and assuming r s = 159.71 Mpc. The median A_D^4 and 
AH values are significantly different from when approximating 
the error on the median as the scatter predicted by the quantilcs 
divided by \/160. Such a discrepancy may be due to our median 
redshift not being exactly z = 0.35 as assumed. It could also be 
due to the breakdown of Taylor assumptions made in deriving our 
fitting model; in the Q m = 0.4 cosmology, our measured a and 
e values deviate substantially from 1 and as shown in the last 
row of Table [2] We again emphasize that while this difference is 
detectable in the median a and e of the mocks, it is not significant 
in each individual mock. 



6 DR7 RESULTS 

6.1 Anisotropic Results 

Now that we have verified the robustness of our techniques 
and obtained a better understanding of the anisotropic sig- 
nal from our mocks, we can proceed to the actual SDSS 
DR7 LRG data. To calculate our fitting model for the data, 
we use the flat ACDM cosmology predicted by WMAP7: 
H = 70.2±1.4 km/s/Mpc, Q. b h 2 = 0.02255±0.054, 9, c h 2 = 
0.1126 ± 0.0036 n 3 = 0.968 ± 0.012 and <t 8 = 0.816 ± 0.024 
l|Komatsu et alJbfJllf ). For the covariance matrix, we again 
use the modified Gaussian covariance matrix discussed in 
H3.3I with the modification parameters derived from the 
mocks and the WMAP7 cosmology. 

The results of our fits are shown in Figure [T5J The 
pre-reconstruction results are in the top row and the post- 
reconstruction results are in the bottom row. The acoustic 
peak appears much sharper after reconstruction, again in- 
dicating the effectiveness of our technique in undoing non- 
linear evolution. This is reflected in the decrease in error on 
a and e by factors of 1.8 and 1.3 respectively after recon- 
struction. The quadrupole near 100/i _1 Mpc scales is much 
closer to after reconstruction which implies that our par- 
tial removal of the Kaiser effect was also successful. The 
deviation from at larger scales again indicates that recon- 
struction is introducing some additional anisotropy. This is 
as expected based on our analysis of the mock catalogues in 
Sj5] (see Figure [6]). 

The first column of Figure [TBI shows the Ax 2 (cv,e) = 
X 2 (q, e) — Xmin distribution measured at various grid points 



in a and e. As described in £13.31 our a grid points are sep- 
arated by 0.0025 in the range 0.7 < a < 1.3 and our e grid 
points are separated by 0.005 in the range —0.3 < e < 0.3. 
The 1 through 6<r confidence levels for a 2D distribution are 
overplotted. The pre-reconstruction results are shown in the 
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Figure 14. Histograms of (a — {a))/cr a (top) and (e — (e))/cr e 
(bottom) after reconstruction. These are a measure of the signal- 
to-noise of our measured a and e values. The best-fit Gaussians 
are overplotted as black lines. We perform a K-S test to see how 
likely these distributions are drawn from a normal distribution. 
The p-values or probabilities are indicated on the plots and im- 
ply that a and e both have finite chances of being drawn from 
Gaussian distributions. This verifies that the standard deviations 
a a and a t we calculate from x 2 ( a > e ) characterize the errors on a 
and e reasonably well. A similar conclusion holds for our a and e 
values before reconstruction. 



top row and the post-reconstruction results are shown in 
the bottom row. As we go to smaller a, the acoustic peak 
in the model is being pushed out to larger scales where the 
errorbars are larger. Hence it is much easier for the fitter to 
hide the peak within the errors. Although we have applied 
a 0.15 prior in log(a) to suppress this unphysical downturn 
in x 2 , the distribution still plateaus in this region. One can 
also see that after reconstruction, the x 2 ( a i e ) distribution 
is much tighter at the center, indicating that the best-fit 
values are much better measured. This corresponds to the 
smaller errorbars we see in a and e after reconstruction. 

The second and third columns of Figure [16] show the 
a and e probability distributions derived from the x 2 grid. 
One can see that both of these are fairly Gaussian so we 
can quantify the errors on a and e as the second moments 
of these distributions, o a and o e . These values are summa- 
rized in Table[3]for both before (top row) and after (bottom 
row) reconstruction. The smaller standard deviations after 
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Figure 15. DR7 fit results before (top row) and after (bottom row) reconstruction. These values imply a 3.6% measurement of Da{z) 
and an 8.4% measurement of H (z) after reconstruction. We see that the acoustic peak has sharpened up significantly after reconstruction 
as expected. The error on a decreases by a factor of 1.8 and the error on e decreases by a factor of 1.3 as a result. The quadrupole is 
nearly at ~ 100/i _1 Mpc after reconstruction, indicating the effectiveness of our Kaiser correction. The deviation from at larger r is 
likely some anisotropy introduced by reconstruction. 



reconstruction accompany the sharpening up of the acoustic 
peak. This corresponds to the tightening of the contours in 
the x 2 distribution shown in column 1. 

Our measured C at , a a and <r e imply correlation coeffi- 
cients of p ae = 0.27 and 0.34 before and after reconstruction. 
These values are slightly larger than the expected p aE ~ 0.21 
from Fisher matrix arguments. However, given the large rms 
of p at from the mocks of ~ 0.35 both before and after re- 
construction, our DR7 results are not significantly different 
from the Fisher matrix prediction. 

Using our measured a, a a , e, a t and C Qe values from 
the DR7 data, we can use Equations ©, ©, (HJ) & © to 
determine Da(z = 0.35) and H(z = 0.35). These results are 
again summarized in Table fn addition, Figure [T71 shows 
the post-reconstruction A^ 2 {Da, H) = \ 2 {Da, H) — Xmin 
contour plot obtained from the a-e grid in Figure 1161 The 
values listed in the table and plotted in the figure were ob- 
tained assuming r a = 152.76 Mpc and using the fiducial 
WMAP7 values: D A j(z = 0M)/r s j = 6.69 and H f {z = 
0.35)r aj / = 12689 km/s. We see that post-reconstruction we 
have a ~ 3.6% measurement of Da(z) and ~ 8.4% measure- 
ment of H(z) from SDSS DR7. Note that our measures of 



Da and H are correlated; this is also clearly evident from 
the contour plot. Before reconstruction, pd a h = 0.30 and 
after reconstruction Pd a h = 0.58. These a r e sim ilar to the 
Pd a h ~ 0.4 predicted bv lSeo fc Eisensteinl (|2007l ). 

We see that our DR7 o a and cr e measurements fall nicely 
within the locus of mock points as shown in Figure[l8]for the 
post-reconstruction case. Note that the mock results shown 
in this figure are identical to Figure[5] The cr e /(l+e)-to-cr a /a 
ratio we obtain is ~ 1.1 before reconstruction and ~ 1.5 after 
reconstruction which is roughly consistent with the Fisher 
matrix prediction of ~ 1.2. Lastly, we note that our oh/H- 
to-ap A I Da ratio is ~ 2 , consistent with the predictions of 
ISeo fc Eisensteinl l|2007h and the assumption that went into 
our Fisher matrix predictions (see Appendix IX)) . 

We again test the robustness of our a and e measure- 
ments to our fitting model. The results are listed in Table [4] 
We see that our a and e measurements are always consistent 
with the results obtained using the fiducial fitting parame- 
ters. Our a a and cr £ measurements show ~ 10% variations 
which are consistent with the differences in scatter between 
the various cases seen in the mocks. 

We also test the robustness of our reconstruction 
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Figure 16. The Ax 2 (ct,e) distribution (column 1) and the derived p(a) and p{e) distributions (columns 2 & 3) for DR7. The pre- 
reconstruction results are in the top row and the post-reconstruction results are in the bottom row. Contour levels corresponding to 1-6<t 
for a 2D distribution are ovcrplottcd. We apply a 0.15 prior in log(a) to suppress the unphysical downturn at low a which corresponds 
to the acoustic peak being pushed out to large r. The errorbars are much larger here and the fitter has an easier time hiding the peak 
inside the errors. The plateauing of the distribution at small a is a result of this. We see that after reconstruction, the contours become 
much tighter. This corresponds to the tightening of p(a) and p(e) after reconstruction seen in columns 2 & 3 due to the sharpening up of 
the acoustic peak. We also see that p(a) and p(e) are nearly Gaussian and hence the second moments a a and <r E characterize the errors 
on a and e well. 



Table 3. Summary of key measurements from DR7 data. Columns 2 and 3 list the a and e values we measure. Column 4 lists the 
covariance between a and e while column 5 lists their correlation coefficient. Columns 6 and 7 list the distance constraints we obtain to 
z = 0.35 from our measured a and e values. Columns 8 and 9 translate these relative distance measures into more tangible quantities 
assuming r s = 152.76 Mpc as in the WMAP7 cosmology. The pre-reconstruction results are listed in the top row (Before) and the 
post-reconstruction results are listed in the bottom row (After). 





a 


e 




Pae 


D A (z)/r s 


H(z)r a 


D A {z) 


H(z) 


Pd a h 














(km/s) 


(Mpc) 


(km/s/Mpc) 




Before 


1.015 ± 0.044 


0.007 ± 0.046 


0.00054 


0.27 


6.751 ± 0.363 


12339 ± 1367 


1031 ± 55 


80.8 ±9.0 


0.30 


After 


1.012 ± 0.024 


-0.014 ±0.036 


0.00029 


0.34 


6.875 ± 0.249 


12895 ± 1079 


1050 ± 38 


84.4 ± 7.1 


0.58 



technique by varying the input parameters and then re- 
performing our fits. The a and e values we measure from 
these tests are also listed in Table [4] Again we see very con- 
sistent a, e, o a and <j e values between the various cases. This 
indicates that our measurements of the acoustic scale and 
anisotropy are robust against reconstruction parameters. 



6.2 Comparison with Past Works 

Three past pa p ers, IGaztanaga. Cabre fe Hul 
IChuang fe Wand (|201lf ) and [Chuang & Wand 



2013), 



20l2), 



have performed anisotrop ic BAO analyses using the DR7 
LRG sample. Of these, IGaztanaga. Cabre fe Hul (|201ll ) 
found a peak in the clustering along the line-of-sight that 
they interpret as a detection of the acoustic peak, but 



iKazin et~al . (2010b) show that the peak is consistent with 
the expected noise. 

IChuang fc Wand i|201ll ) measure D A (z = 0.35) and 
H(z — 0.35) by fitting the 2D correl ation function of 
the L RGs before reconstruction, whereas IChuang fc Wand" 
(|2012T ) fit the monopole and quadrupole. The former mea- 
sures D A (z = 0.35) = 10481^ Mpc and H(z = 0.35) = 
82.1^4g km/s/Mpc at z = 0.35 while the latter measures 
D A {z = 0.35) = 1057l|? Mpc and H(z = 0.35) = 79-6±|;f 
km/s/Mpc. These values are all consistent with our pre- 
reconstruction measurements, however, the magnitudes of 
the errors are slightly different. Since our treatments of the 
covariance matrix differ, this is not surprising. 

A significant differe nce between our analy s is and that of 
IChuang fc Wand (|201ll ) and lChuang fc Wand (|2012T ) is that 
they use a Markov Chain Monte Carlo approach over the pa- 
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Table 4. DR7 fitting results for various models. The model is given in column 1. The measured a values are given in column 2 and the 
measured e values are given in column 3. The x 2 /dof is given in column 4. 



Model 


a 


e 


X 2 /dof 


Redshift Space without Reconstruction 


Fiducial [/] 


1.015 ±0.044 


0.007 ± 0.046 


89.60/90 


Fit w/ (E ± , E||) -> (8, 8)/i _1 Mpc. 


1.012 ±0.045 


0.009 ± 0.044 


89.77/90 


Fit w/ E s -> Oft _1 Mpc. 


1.018 ±0.040 


0.007 ± 0.040 


89.60/90 


Fit w/ A 2 (r) = poij/2. 


1.018 ±0.043 


0.013 ± 0.046 


91.42/91 


Fit w/ A 2 (r) = polyi. 


1.015 ±0.044 


0.006 ± 0.047 


89.58/89 


Fit w/ 30 < r < 200ft _ 1 Mpc range. 


1.018 ±0.039 


0.004 ± 0.043 


105.03/104 


Fit w/ 70 < r < 200ft - 1 Mpc range. 


1.016 ±0.050 


0.008 ± 0.050 


82.43/76 


Fit w/ 50 < r < 150ft" 1 Mpc range. 


1.019 ±0.042 


0.001 ± 0.049 


47.10/58 



Redshift Space with Reconstruction 



Fiducial [/] 




1.012 ±0.024 


-0.014 ±0 


036 


62.53/90 


Fit w/ 


(E X ,E||) -> (2,4)ft- 1 Mpc. 


1.012 ±0.025 


-0.014 ±0 


036 


62.48/90 


Fit w/ 


E s -> 


Jft^Mpc. 


1.013 ±0.021 


-0.013 ±0 


029 


61.83/90 


Fit w/ 


A 2 (r) 


= poly2. 


1.013 ±0.025 


-0.011 ± 


036 


65.61/91 


Fit w/ 


A 2 (r) 


= polyA. 


1.013 ±0.025 


-0.011 ± 


036 


61.92/89 


Fit w/ 


30 < r 


< 200ft~ 1 Mpc range. 


1.014 ± 0.023 


-0.013 ±0 


033 


68.39/104 


Fit w/ 


70 < r 


< 200ft- 1 Mpc range. 


1.012 ±0.027 


-0.016 ±0 


040 


54.50/76 


Fit w/ 


50 < r 


< 150ft - 1 Mpc range. 


1.017 ± 0.023 


-0.009 ±0 


034 


31.95/58 


Recon. 


w//9- 


-> 0.24. 


1.014 ±0.024 


-0.016 ±0 


035 


68.77/90 


Rccon. 


w//3- 


-> 0.36. 


1.013 ±0.024 


-0.013 ±0 


036 


67.05/90 


Recon. 


w/ 6 - 


> 1.8. 


1.014 ± 0.025 


-0.017 ± 


036 


66.75/90 


Recon. 


w/ b - 


> 2.6. 


1.015 ±0.024 


-0.012 ±0 


035 


77.09/90 


Recon. 


w/ Wiener Filter. 


1.012 ±0.025 


-0.014 ±0 


035 


61.23/90 
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Figure 17. The post-reconstruction SDSS DR7 Ax 2 {D A ,H) = 
X 2 {Da, H) — Xmin distribution with 1 and 2a contours overplot- 
ted. We measure D A (z) = 1050 ± 38 Mpc and H(z) = 84.4 ±7.1 
km/s/Mpc at z = 0.35. The tilted elliptical contours clearly in- 
dicate a correlation between our D A and H measurements. The 
correlation coefficient is Pd a h = 0.58. 



rameter space {Da(z), H(z), /3, Q.mh 2 , Q b h 2 , n 3 , E s , E ni } to 
derive their Da and H measurements. The inclusion of /? 
and E s in their parameter estimation without any other 
marginalization (such as A(r)) to allow for inaccuracies in 
their model indicates that their method also attempts to ex- 
tract anisotropic information in the broadband shape of the 
correlation function. This differs from the method we present 
whic h only utilizes the an isotropy of the BAO signal. How- 
ever. iReid fc White! |201ll ) has shown that accurately mod- 
eling broadband redshift-space distortions is a challenging 
theory problem and requires going beyond simple E s and /3 
models. They demonstrate that neglecting bispectrum and 
higher order terms from the real to redshift space transfor- 
mation results in models that are not accurate enough at 
our current le vels of observational pre cision. Hence, as also 
pointed out in lChuang fc Wand (|2012l ). using these simple, 
inexact models without additional marginalization can re- 
sult in biased measurements of Da and H. This may also 
lead to the difference in their measured errors compared to 
ours. 

Our DR7 q measurements are con s istent with the 
monopole-only measurements of IXu et al.l {2012), however, 
our errorbars on a are a factor of 1.25 larger both before and 
after reconstruction. Although this is a very small change in 
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Figure 18. Post-reconstruction a a versus cr e for the mocks with 
the DR7 point overplotted as the black star. Note that the mock 
points arc identical to Figure [9] We see that our DR7 mea- 
surement falls nicely within the locus of mock points. The DR7 
cr e /(l + t)-to-a a /a ratio is ~ 1.1 before reconstruction and ~ 1.5 
after reconstruction, roughly consistent with the Fisher matrix 
prediction of ~ 1.2. 



absolute terms, it is still worthy of some investigation. The 
top left panel of Figure [19] shows the a a values we measure 
for our full monopole+quadrupole fits as described in £13.31 
from the post- reconstru c tion m ocks versus the monopole- 
only results of IXu et al.l jjoij ). The DR7 results are over- 
plotted as the black star. We see that, in general, the mocks 
have larger a a values in our full monopole+quadrupole fits. 
The scatter at large a a where the acoustic scale is not well 
measured is also significantly bigger. Our DR7 a a measure- 
ment lies within the locus of mock points and hence the 
increase we see is consistent with the mocks. 

Our monopole+quadrupole full fits h ave severa l differ - 
ences relative to the monopole-only fits of IXu et al.l (|2012l ). 
First, the covariance matrix is expanded to include the 
quadrupole-quadrupole and monopole-quadrupole covari- 
ances. Second, we introduce /3 as a fitting parameter. Third, 
we includ e FoG (i.e. E s = 4h~ Mpc) in our fitting model; 
IXu et all (|2012l ) have E s implicitly set to 0/i _1 Mpc. Since 
E s can induce some smearing of the BAO, the E n i value 
we use in the full fits is correspondingly smaller (3/i _1 Mpc 
versus 4/j~ 1 Mpc in the monopole-only case). In addition, 
in our pre-reconstruction fitting model, we introduce a 
non- isotropic E n i. To understand which of these steps in- 
duces the greatest change in a a , we start by fitting the 
monopole+quadrupole using the new covariance matrix and 
gradually add in the other changes. Before we describe our 
results, we again stress that the changes in a a we see are 
very small and require probing some subtleties in our mod- 
els to understand. Our measurements of a and a a are still 
reasonably robust against fitting parameters in a single DR7 
realization as shown in Tables [2] & [4] 

The median changes in a a (relative to the monopole- 
only case) as we add in more elements of the fitting are 
listed in Table[S]for the pre- and post-reconstruction mocks. 
The total change between the monopole-only fits and the 
full monopole+quadrupole fits is listed in the last row of 
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Figure 19. The a a values measured from the post-reconstruction 
mocks for various fitting models. The DR7 results are over- 
plotted as the black stars and fall within the locus of mock 
points, (top left) The a a values measured through the full 
monopole+quad rupole fits versus the monopole-only results of 
IXu et aU |2pi2). One can see that the full fits have a a that 
are larger on average. There is also considerable scatter at 
large a a where the acoustic scale is not as well measured. We 
emphasize however, that these variations in a a are incredibly 
small and do not significantly affect our results, (bottom left) 
a a from monopole+quadrupole fits with S s = 4fc _1 Mpc and 
S n i = 3/i -1 Mpc versus those from the monopole-only fits with 
E s = 0h _1 Mpc and S nl = 4ft, _1 Mpc. The degradation in (J a 
is obvious and accounts for half the increase in a a relative to 
the monopole-only case. This again suggests that there is some 
mismatch between our fitting model and the data, (top right) cr a 
from full monopole+quadrupole fits including e versus those from 
the monopole+quadrupole fits with e = 0. The introduction of e 
appears to cause most of the scatter in a a and is responsible for 
the other half of the <r Q increase relative to the monopole-only 
case, e is known to have some correlation with a so this is not 
surprising. 

this table. We find that changing the covariance matrix (first 
row) increases a a by a small amount and adding f3 (second 
row) does not further degrade the errors. Introducing E s 
and the accompanying change in E n i (third row) appears 
to be a major contributor to the degradation of a a . This 
increases the median a a by about half the total. Combining 
the (3 fitting and the changes in E s and E n i (fourth row) 
shows little additional degradation above the previous case. 
Note that this corresponds to using the fiducial model with 
e fixed at 0. Finally, as mentioned above, the last row adds 
in e fitting and corresponds to the fiducial model. We see 
that this step causes the other half of the total increase in 
a a . It also introduces a significant amount of scatter in a a . 

The steps that contribute the most to the a a in- 
crease are shown in the bottom left and top right pan- 
els of Figure [19] for the post-reconstruction case. In the 
bottom left we have plotted the a a values measured from 
monopole+quadrupole fits with E s = 4ft" 1 Mpc and E n i = 
versus those measured from the monopole-only 
fits with E s = 0/i _1 Mpc and E n i = 4/i _1 Mpc. The offset 
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Table 5. Changes in a a relative to the monopole-only fits of IXu et al,l l|2012h that arise when we introduce new fitting elements. In 
the first row we introduce the quadrupole fitting with the new combined monopole+quadrupolc covariancc matrix. The second row 
introduces our fitting in addition to the new covariance matrix while the third row introduces changes in S s and S n i instead. The 
fourth row combines the previous two and corresponds to fitting with the fiducial model while forcing e = 0. The last row introduces e 
fitting and corresponds to our fiducial model results. 



Parameters Before Recon. After Recon. 

Acr Q Qtls Acr a Qtls 

New dj 

New C\ j , adding f) fit 
New dj, S s = 4/i _1 Mpc, new S nl 

New C'ij, adding /3 fit, S s = Ah~ 1 Mpc, new S n i (fiducial model w/ e = 0) 
New Cij, adding /3 & e fits, £„ = 4/i -1 Mpc, new S n [ (fiducial model) 



0.0007 
0.0007 
0.0015 
0.0018 
0.0036 



+0.0014 
-0.0015 
+0.0030 
-0.0016 
+0.0017 
-0.0016 
+0.0024 
-0.0017 
+0.0040 
-0.0044 



0.0005 
0.0006 
0.0013 
0.0014 
0.0022 



+0.0020 
-0.0016 
+0.0018 
-0.0019 
+0.0021 
-0.0016 
+0.0014 
-0.0016 
+ 0.0039 
-0.0027 



between the two is obvious and again suggests that our FoG 
model is not perfectly matched to the data. 

The top right panel shows the a a values measured from 
the mocks through the full monopole+quadrupole fits versus 
the monopole+quadrupole fits with e fixed at 0. In addition 
to the obvious offset, we also see the appearance of signif- 
icant scatter, e has small but non-zero correlation with a, 
implying a slight degeneracy between these 2 parameters. It 
is not surprising that this extra covariance may increase a a 
and its scatter. Again, this small degradation is not of great 
concern at our current levels of statistical precision. 



7 COSMOLOGICAL IMPLICATIONS 

In this section we will place our measurement of e within 
the context of current cosmological constraints. To build 
more intuition for how a(z) and e(z) vary as we change the 
amount of curvature or the nature of dark energy, we look 
to Figure [20l The left panel shows a as a function of redshift 
for a cosmology that has positive curvature (Ok = 0.1), dark 
energy that is not a cosmological constant (wo = —0.9) and 
time- varying dark energy (w a = 0.5). The analogous plot for 
e is shown in the right panel. Here we have again taken the 
fiducial cosmology to be the flat ACDM cosmology predicted 
by WMAP7. As we vary the cosmology, we fix Q m h 2 and 
the distance to the last scattering surface (i.e. the distance 
to z = 1089, the redshift of recombination). This guarantees 
that the sound horizon and the CMB remain approximately 
unchanged in all the plotted cosmologies. 

Recall that if the fiducial cosmology matches the true 
cosmology of the universe, then we would expect a = 1 and 
e = 0. We see that introducing curvature and altering the 
nature of dark energy both perturb a away from 1 the most 
at low redshift. However, at higher redshift, it becomes in- 
creasingly difficult to distinguish between non-cosmological 
constant models and time- varying dark energy models using 
measurements of a. The opposite is true for e. We see that 
the effects of adding curvature or changing the properties of 
dark energy are most prominent at larger redshifts, peak- 
ing at z ~ 1. This suggests that to exploit the anisotropic 
BAO signal, we gain more leverage by going to higher z. 
However, we also see that even the maximum difference in e 
between the w a — 0.5 and SIk = 0.1 cosmologies is smaller 
than our current error on e which indicates that we are not 
able to distinguish between these cosmologies using our DR7 
measurement. 



Figure [5T] shows our DR7 measurement of e overplot- 
ted on constraints derived from cosmological Markov Chain 
Monte Carlo (MCMC) sampling. The MCMC method com- 
putes the likelihood that a set of input cosmological param- 
eters fits distance measures from Cosmic Microwave Back- 
ground (CMB) observations at high redshift, and Type la 
supernova (SN) and BAO observations at low redshift. The 
number of steps in the chain spent exploring a certain region 
in the cosmological parameter space is proportional to the 
likelihood of that region representing the true cosmology. 
Hence, we can infer Da(z) and H(z) at each step in the 
chain to compute e(z) relative to some fiducial cosmology 
(WMAP7 in our case). At each z, we can measure the mean 
and rms of the e distribution which can then be compared 
to our DR7 measurement at z = 0.35. 

The grey regions in Figure |2"T1 are derived from MCMC 
chains exploring a cosmology which is allowed to have cur- 
vature and varying dark energy with equation of state 
w(a) = wo + (1— a)w a where a is the scale factor. Note that 
this is the most generalized and least-constraining cosmol- 
ogy that is typically tested and is the reference cosmology 
defi ned by the Dark Ener gy Task Force (the DETF cosmol- 
ogy; [Xlbrechtj[riil] [2006|). These chains we re computed us- 
ing C osmoMC, a standard MCMC sampler (jLewis fc Bridle! 
|2002| ). The dark grey and light grey regions correspond to 
the la limits calculated from chains using CMB+BAO data 
and CMB+BAO+S N data respectively . The CMB data i s 
taken from WMAP7 (|jarosik et alj201ltlLarson et al.ll201ll ) 
and the SN distance contraint s are taken from th e Super- 
nova Legacy Survey 3 (SNLS3; IConlev et alJlioTH ). 

The BAO distance measure used in the chains is 
based on our DR7 anisotropic measurements, however, we 
marginalize over e and only use the remaining isotropic 
or spherically-averaged information encoded in a. This is 
equivalent to using p(a), the probability distribution for the 
spherically-averaged distance Dy to 2 = 0.35 (plotted in 
Figure [T6)l . obtained by applying Equation (|53[) to our mea- 
sured p(a, e). Figure |2"T1 overplots our DR7 e measurement 
on the e constraints obtained through our chains. One can 
see that our DR7 measurement of e is consistent with the 
CMB+BAO and CMB+BAO+SN constraints at z = 0.35. 
In the CMB+BAO case, our DR7 e errorbar falls within 
the e constraints obtained by using the spherically-averaged 
distance measure. This suggests that the DR7 anisotropic 
information we measure offers additional constraints on the 
allowed cosmological parameter space. 
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Figure 20. The expected variation in a(z) (left) and t(z) (right) as we open up curvature or allow non-cosmological constant or time- 
varying dark energy. We have taken the fiducial cosmology to be the flat, ACDM cosmology predicted by WMAP7 as usual. One can 
see that curvature and dark energy properties affect a more at low z. At high z it becomes increasingly difficult to distinguish between 
non-cosmological constant and time- varying dark energy models using measurements of a. However, we see that t is affected by curvature 
and dark energy the most at higher z, peaking at z ~ 1. This suggests that the anisotropic BAO signal is stronger and therefore offers 
more constraining power at higher rcdshifts. 



Table 6. Cosmological parame ters measured from CosmoMC chains. The first line listed for each cosmology is the DR7 monopole-only 
result from lMehta et alj l|2012l ). The second (bold) line is derived from our DR7 a and t measurements. 



Cosmology 


Datasets 


£^777, h 


£~^m 


Ho 




w 


W a 


oCDM 


CMB+BAO 


0.1333(53) 
0.1342(50) 


0.278(15) 
0.285(16) 


69.3(16) 
68.7(18) 


-0.004(5) 
-0.003(5) 






wCDM 


CMB+BAO 


0.1349(57) 
0.1347(57) 


0.285(25) 
0.289(26) 


69.0(39) 
68.4(38) 




-0.97(17) 
-0.95(17) 




OWQWaCDM 


CMB+BAO+SN 


0.1346(53) 
0.1352(52) 


0.276(15) 
0.280(17) 


69.9(19) 
69.6(21) 


-0.010(7) 
-0.010(7) 


-0.90(16) 
-0.90(16) 


-1.30(99) 
-1.32(100) 



The bold text in Table [6] lists results from CosmoMC 
chains run using our DR7 a and e measurements. In the 
table, oCDM refers to a cosmology in which Ok is al- 
lowed to vary (i.e. the universe is allowed to be curved), 
wCDM refers to a cosmology in which dark energy is al- 
lowed to vary from a cosmological constant and owow a CDM 
is the DETF cosmology. Here we have combined our full 
anisotropic BAO measurements at z = 0.35 with ob serva- 
tions of the CMB from WMAP7 (|Komatsu et al.ll201ll ) and, 
in the c ase of the DETF cosmology, supernovae from the 
SNLS3 i|Conlev et"al"1l201ll). W e compare these with the re- 
sults from Mehta et alj ( 20121 ) (non-bold text in the table) 
and find consistent measurements and errors for the cos- 
mologies tested. In addition, using the e-marginalized p(a) 
distribution obtained by applying Equation ([53} to p(a, e) 
gives similarly consistent values of the cosmological param- 
eters. However, as expected, the errors can be ~ 10 — 20% 
higher in this case since we are not including the DR7 
anisotropic constraints. 

Recall that the errors on a we mea s ure fr om DR7 are 
larger than those m easured bv lXu et al.1 (|2012h and used in 
iMehta et al](|2012D . The fact that our chains produce similar 
errors on the cosmological parameters implies that we are 
obtaining some constraint from e. Since the fully anisotropic 
analysis presented here is a more careful treatment of BAO 
measurements, the similarity in results confirms the validity 



of past monopole-only BAO analyses. We also note that the 
redshift of DR7 (z — 0.35) is relatively low and as shown 
in Figure 1201 the main constraining power of e comes in 
at larger redshifts. This implies that anisotropic analyses 
should be more powerful than monopole-only analyses as 
we go to higher z. 

The future is bright w ith the SDSS DR9 CMASS sam- 
ple |Anderson et al.ll2012i ) now in hand. The galaxies in this 
dataset are denser and more abundant than DR7, and also 
at higher redshift (z = 0.57). Therefore, the DR9 CMASS 
anisotropic BAO signal should be less noisy and more promi- 
nent, implying a much tigher constraint on e. 



8 CONCLUSIONS 

The differential clustering along the line-of-sight and trans- 
verse directions that arise from assuming the wrong fidu- 
cial cosmology can be used to directly constrain the angular 
diameter distance Da(z) and the Hubble parameter H(z). 
This anisotropy can be measured from the BAO signal in 
the monopole and quadrupole moments of 2-point statistics 
such as the correlation function studied in this work. 

We have presented measurements of the anisotropic 
BAO signal (e) from the SDSS DR7 LRG sample including 
density-field reconstruction. We measured a = 1.012 ±0.024 
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Figure 21. Our DR7 e measurement at z = 0.35 overplottcd 
on the la regions predicted by MCMC chains computed via com- 
bining CMB+BAO (dark grey) and CMB+BAO+SN (light grey) 
datasets. The solid black line corresponds to the mean e(z) from 
the CMB+BAO chain and the dashed black line is the analogue 
for the CMB+BAO+SN chain. Here we have assumed a WMAP7 
fiduc ial cosmology. The p lotted chains assumed a DETF cosmol- 
ogy l|Albrecht et al]|200rj ) in which the universe is allowed to be 
curved and the dark energy equation of state takes on the form 
w(a) = w +(l-a)wg.The CMB data u sed is taken from WMAP7 
jjarosik et al ] |201ll ; lLarson et alj201ll) and the SN distance con- 
straints used are taken from the SNLS3 ijConlev et al 1 l201lh . The 
BAO distance constraints comes from applying Equation < 1 5 3 D to 
the DR7 p(a,e) distribution in order to obtain p(ct), collapsing 
our 2D constraint into a single constraint on the spherically aver- 
aged distance Dy to z = 0.35. This is equivalent to marginalizing 
over e. One can see that our DR7 e measurement overlaps the re- 
gions predicted by our chains very well. In addition, we see that 
our errorbars on e are smaller than the region predicted by using 
the spherically averaged constraint p(a) in the CMB+BAO case. 
Hence the anisotropic information encoded in our measurement 
of e provides additional constraints on the parameter space of 
allowed cosmologies. 



and e = -0.014+0.036 which translate into D A (z = 0.35) = 
1050 + 38 Mpc and H(z = 0.35) = 84.4+7.1 km/s/Mpc as- 
suming r s — 152.76 Mpc. Note that these measurements 
of Da and H are correlated with Pd a h = 0.58. We have 
demonstrated that the methods for extracting e outlined in 
this paper are robust and applicable to future anisotropic 
BAO studies. 

We have given a detailed account of the theoretical 
background motivating the origin of the anisotropic signal 
and a parameter, e, for measuring it. An in-depth look at the 
fitting model and method we use to extract the anisotropic 
signal is also given. We find that our model parameters have 
different morphological structures in their derivatives from 
e, although they can still be partially degenerate with each 
other. These minor degeneracies appear to introduce a small 
bias in e at the 0.2% level, far below our current level of sta- 
tistical precision. 

We apply density field reconstruction and test the ro- 
bustness of our measured a and e against changes in the 
reconstruction parameters using 160 LasDamas mock cat- 
alogues. We find that reconstruction appears to introduce 
some anisotropy into the quadrupole, however this is ad- 
equately accounted for by our A2(r) nuisance parameters. 
We then perform the same robustness checks on our fitting 



model using the mock catalogues. Similar tests were also 
performed on the DR7 data returning consistently robust 
results. We demonstrate that a and e have near-Gaussian 
posteriors. Hence estimating their errors (a a and cr e ) as the 
second moments of their respective probability distributions 
is reasonable. The a a , cr e and p ae values obtained from the 
mocks and the DR7 data are mostly consistent with Fisher 
matrix predictions. 

We find that in the mocks and the DR7 data, our a 
error estimates are slightly larger than those obtained when 
only the monopole is fit. This small increase does not de- 
tract significantly from the overall robustness of our mea- 
surements which we verify as discussed above. About half of 
this increase is a result of including an FoG model in our full 
monopole+quadrupole fits while the other half arises from 
fitting for e. This first point suggests that our FoG model 
does not match the data perfectly and may induce slight bi- 
ases into our measurements. Given the non-zero covariance 
between e and a the second point is not surprising. The be- 
haviour of DR7 falls completely within the locus of mock 
points and is therefore not unusual. 

Our DR7 measurements of Da and H before reconstruc- 
tion a re co nsistent with thos e obta ined bv lChuang fc Wand 
(|201lh and lChuang fc Wand (|2012h using the same dataset. 
The errors we measure, however, differ slightly due to differ- 
ences in our analysis techniques. The errors on a and e we 
measure from the DR7 data are consistent with the scatter 
from the mocks. In addition, our cr e /(l + e)-to-a a / a ratio 
agrees reasonably well with Fisher matrix predictions and 
our oh / H-to-ap a /Da ratio is ~ 2 and i n goo d agreement 
with the predictions of ISeo fc Eisensteinl (|2007l ). 

Our post-reconstruction DR7 e measurement agrees 
well with the predictions from current datasets. Compar- 
ing the cosmological parameters we obtain using our a and 
e measurements w ith the monopole-only measurements in 
iMehta et af] (|2012Tl yields consistent values and errors. This 
suggests that although our errors on a are larger than the 
monopole-only case, the e measurement is providing some 
additional constraint on the cosmological parameters. The 
validity of previous monopole-only analyses is also confirmed 
by the similarity between those cosmological results and the 
ones obtained here through a more careful, fully-anisotropic 
analysis of the BAO. 

We find that the anisotropic signal is stronger at higher 
redshifts which suggests that its constraining power will be- 
come more apparent at higher z. The recently obtained 
SDSS DR9 CMASS dataset has a higher number density 
than the DR7 LRG sample, contains more galaxies and is at 
higher redshift (z = 0.57). The basic theory and methodol- 
ogy presented in this work should serve as a foundation for 
obtaining a much better detection of e, and subsequently 
D A , H and other cosmological parameters from CMASS. 
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APPENDIX A: FISHER MATRIX 
PREDICTIONS 

Using a Fisher matrix formalism, it is possible to derive the- 
oretical predictions for and expected correlations between 
the variances and covariances of a, e, Da and H. The val- 
ues derived here are utilized in the text as sanity checks for 
the values we measure. 

We begin with the matrix equation 



is typically ~ 2 (|Seo fc Eisensteir] 120071 ') . The correlation 
coefficient Pd a h is predicted to be ~ 0.4. Hence, we have 



— °"log a — 0.73<Tl O g o A 

a 

= Clog(l+e) = 0.86(Ti og £) A , 



1 + e 

which implies the ratio 

to ~ 1.2. 

1 + e a 

The correlation coefficient between a and e is 
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Using Equations (|A11[) & (|A12[) and plugging in the assumed 
values of / and Pd a h gives 



0.21. 
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Note that this is essentially the inverse process to Equation 
0. Expanding we get 
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Plugging in the relevant derivatives from Equations Q & 
(0} we get 
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where af ogy = 3f. 

The correlation coefficient between Da and H is 
Pd a h = &d a h I & d a ®h - If we write / = a log H / m og d A , then 
we have 



(1 + e)* 
a(l + e) 
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Note that / is just the ratio of a h / 'H-to-a d a / 'D a which 



